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This  thesis  defined  a  fundamental  approach  for  aircraft  wingbox  design  t^ipropriate  for  an  introductory 
course  in  aircraft  structures  based  upon  material  strength  and  stiffness  requirements.  The  process  developed 
sought  to  encompass  major  conceptual  engineering  desip  considerations  that  ranged  from  load  estimation  at 
various  points  in  the  subsonic  flight  envelope,  to  initial  structural  sizing  and  layout.  The  goal  was  to  present  a 
process  that  could  be  readily  conducted  via  hand  calculations  and  applied  by  any  student  entering  basic  aircraft 
structures  design. 

The  sequence  of  analysis  began  with  application  of  a  comprehensive  panel  code  method  developed  by 
NASA  Ames  Research  Center  known  as  PMARC.  Loads  obtained  from  the  code  were  then  used  to  foimulate 
a  strength  of  materials  study  of  the  structure  subjected  to  combined  bending,  shear  and  torsion.  The  static  load 
approach  allowed  initial  estimation  of  component  sizing  based  upon  materi^  or  buckling  allowable  stress 
selection.  Finally,  the  study  demonstrated  a  strength  to  weight  ratio  comparison. 

Several  calculation  examples  and  computer-based  spreadsheets  were  prepared  for  rapid  analysis  of  multiple 
option  desip  scenarios.  Since  the  study  employed  an^ysis  methods  that  could  be  performed  without  the  aid 
of  a  finite  element  routine  or  extensive  computer  programming  knowledge,  it  serves  as  a  good  introduction  for 
the  entry  and  intermediate  level  structural  engineer. 
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ABSTRACT 


This  thesis  defined  a  fundamental  approach  for  aircraft  wingbox  design  appropriate 
for  an  introductory  course  in  aii  craft  structures  based  upon  material  strength  and  stifihess 
requirements.  The  process  developed  sought  to  encompass  miyor  conceptual  engineering 
deagn  conaderations  that  ranged  from  load  estimation  at  various  points  in  the  subsonic 
flight  envelo^  j  u  al  structural  sizing  and  layout.  The  goal  was  to  present  a  process 
that  could  be  readily  conducted  via  h&nd  calculations  and  applied  by  any  student  entering 
basic  aircraft  structures  design. 

The  sequence  of  analyis  began  with  application  of  a  comprehensive  panel  code 
method  developed  by  NASA  Ames  Research  Center  krown  as  PMARC.  Loads  obtained 
fi-om  the  code  were  then  used  to  formulate  a  strength  of  materials  study  of  the  structure 
subjected  to  combined  bending,  shear  and  torsion.  The  static  load  approach  allowed  initial 
estimation  of  component  sizing  based  upon  material  or  buckling  allowable  stress  selection. 
Finally,  the  study  demonstrated  a  strength  to  weight  ratio  comparison. 

Several  calculation  examples  and  computer-based  spreadsheets  were  prepared  for 
rapid  analysis  of  multiple  option  design  scenarios.  Since  the  study  employed  analysis 
methods  that  could  be  performed  without  the  aid  of  a  finite  element  routine  or  extensive 
computer  programming  knowledge,  it  serves  as  a  good  introduction  for  the  entry  and 
intermediate  level  structural  engineer. 
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I.  INTRODUCTION 


This  thesis  provides  a  systematic  approach  to  the  preliminary  design  of  an  aircraft 
v^ngbox  structure  that  employs  a  two  spar  support  base.  The  method  chosen  was  based 
upon  an  applied  elasticity  approach  as  described  by  Allen  and  Haisler  [Ref  1  :pp.  145-233] 
using  strength  and  stiflhiess  specifications  as  the  design  target.  The  intent  of  this  efibrt 
was  to  provide  insight  into  the  flow  of  engineering  design  decisions  that  begin  in  load 
determiiuition  and  end  with  preliminary  structural  layout.  As  a  result  of  this  work,  several 
study  sections  were  developed  for  incorporation  into  the  current  aeronautical  engineering 
curriculum  at  the  Naval  Postgraduate  School  (NPS)  to  be  presented  in  the  undergraduate 
and  graduate  core  structural  engineering  courses. 

The  organization  of  this  thesis  begins  with  determination  of  the  static  load  boundary 
condition  applied  to  an  airfoil  of  arbitrary  geometry.  To  solve  the  low  qjeed  problem  for 
designs  encountering  flight  speeds  up  to  and  including  Mach  number  0.6,  a  potential  flow 
panel  code  method  was  chosen  that  allowed  user  definition  of  wing  geometry,  and  total 
aircraft,  if  desired.  The  code  was  used  in  this  research  to  calculate  spanwise  force  and 
moment  information;  however,  it  is  also  capable  of  providing  pressure  distribution 
information  along  the  surface  of  the  structure  as  well. 

Once  the  loads  were  obtained  along  the  span  from  the  code,  they  were  used  as  a 
series  of  discrete  point  loads.  From  this  information  traditional  bending  and  shear 
diagrams  were  constructed  for  use  in  the  strength  analysis.  As  part  of  the  instructional 
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module  for  the  Introductory  Aircraft  Structures,  AA  2021,  a  teaching  guide  and 
automated  spreadsheet  routine  were  developed  to  aid  in  presenting  tliis  material  in  a  four 
hour  lecture  series. 

The  thesis  continues  with  a  review  of  the  buckling  stability  of  a  wing  cover  subjected 
to  the  stresses  projected  from  a  shear  flow  solution.  Using  the  methods  described  in 
Sechler  and  Dunn,  an  analysis  of  a  stiffened  wing  section  was  conducted  to  determine 
potential  stringer  size,  geometry  and  total  number  required  to  support  the  structure  at  a 
given  cross-section  [Ref  2;pp.  1 54-258].  As  these  chapters  are  presented,  several  design 
questions  arise  from  the  structural  engineer's  perspective  that  are  presented  and  answered. 

At  the  conclusion  of  the  research,  several  follow-on  topics  were  presented  for  future 
thesis  work  by  students  at  NFS.  Since  this  woilc  does  support  an  ongoing  design  effort 
that  is  a  fundamental  segment  of  the  curriculum  for  all  students  within  the  department,  this 
work  can  serve  as  a  valuable  reference  for  use  in  future  aircraft  design  reviews  and 
competitions. 
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11.  LOAD  ANALYSIS 


A.  REQUIREMENTS 

The  very  nature  of  design  requires  tools  that  support  creative  and  interactive 
development  of  aircraft  structures.  Interactive  flexibility  demands  that  the  de«gner  must 
be  able  to  quickly  analyze  several  options  feeding  back  to  arrive  at  a  configuration  that 
best  suits  the  requirements  at  hand.  The  boundary  condition  imposed  upon  an  aircraft 
structure  is  the  result  of  an  external  pressure  distribution,  internal  structural  weight  and 
internal  pressure  distribution  associated  with  a  given  flight  envelope.  These  components 
of  the  overall  boundary  condition  can  then  be  idealized  in  a  number  of  ways.  Any  method 
chosen  to  idealize  the  boundary  condition  to  allow  solution  of  its  impact  on  the  strength  of 
a  structure,  necessitates  the  capability  to  adapt  to  arbitrary  geometries  and  provide 
meaningful,  accurate  data  in  a  reasonable  amount  of  time. 

Several  Computational  Fluid  Dynamics  (CFD)  methods  exist  today  that  provide  very 
reliable  data  for  nearly  all  aspects  of  the  flight  envelope.  Unfortunately  these  methods  also 
require  extensive  training  and  experience  in  grid  generation,  input  data  preparation, 
interpretation  of  results,  etc.;  in  addition  to  tremendous  computer  processing  time. 

Another  drawback  to  CFD  is  the  time  associated  with  development  of  an  individually 
tailored  analysis  grid  which  can  require  several  weeks  of  effort. 

Earlier  estimation  methods  employing  converrion  of  2-d  airfoil  data  into  3-d  wing 
results,  as  described  by  Nicholai  [Ref  3;pp.  2-l:2-9]  and  Raymer  [Ref  4;pp.  342-344] 
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tend  to  leave  many  students  searching  for  specific  solutions  to  problems  that  are  not 
covered  by  generic  tables  or  charts 

B.  REVIEW  OF  AVAILABLE  METHODS 

For  the  reasons  listed  above,  the  CFD  approach  was  rejected  for  use  in  an 

introductory  to  intermediate  level  approach  to  load  determination.  Since  most  tactical 
naval  aircraft  operate  in  both  the  supersonic  and  subsonic  flight  regimes,  methods  of  load 
approximation  are  required  in  both  arenas.  The  problem  of  estimating  supersonic  loads 
was  believed  to  be  reasonably  approximated  by  the  application  of  2-d  gas  dynamics  to 
obtain  the  external  resultant  pressure.  Since  all  supersonic  aircraft  also  operate  in  the 
subsonic  regime,  the  load  calculation  method  chosen  for  this  work  was  for  those  loads 
associated  with  flight  operations  in  the  0.0-0.6  Mach  number  regime.  One  area  this  thesis 
did  not  address  concerns  the  load  approximation  for  flight  in  the  transonic  regime.  For 
first  order  approximations  at  those  speeds,  the  reader  is  encouraged  to  review  the  methods 
depicted  in  Raymer  [Ref  4:pp  293-297]. 

Several  analytical  and  empirical  approximations  in  the  subsonic  range  have  been 
developed  throughout  the  years  to  aid  the  designer  in  a  quick  approximation  of  the  loads 
experienced  along  the  span  of  a  wing.  Two  of  the  more  recognized  approaches  include 
DATCOM  and  Schrenk's  approximation.  Unfortunately  during  previous  work,  several 
students  found  these  methods  lacking  in  sufficient  detail  when  confironted  with  changing 
geometries,  (consider  the  case  of  high-lift  devices),  or  actual  changes  in  airfoil 
cross-section,  (i.e.,  wash-out).  However  these  methods  did  prove  to  be  useful  in 
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validating  computer  solution  methods.  A  key  design  lesson  learned  in  the  application  of 
any  computer  solution  is  a  thorough  validation  of  expected  results  from  a  second 
approach  at  one  or  more  known  benchmarks. 

Students  are  exposed  to  potential  flow  theory  early  in  the  aeronautical  engineering 
curriculum  at  NPS  and  generally  understand  the  underlying  principles  of  most  pane!  code 
solution  methods.  Prior  to  replacement  of  the  department's  VAX  computer  network 
system,  a  panel  code  produced  by  NASA  Langely  known  as  SUB  had  been  used 
extensively  as  a  design  tool  by  many  students.  During  the  interim  between  the  loss  of  the 
VAX  and  complete  system  integration  of  a  UNIX  based  Silicon  Graphics^  workstation 
network,  the  department  sought  to  incorporate  a  new  panel  code  routine  for  use  by  the 
aerospace  engineering  student  body. 

C.  PMARC 

During  1993  LCDR  Dave  Porter,  USN,  obtained  a  panel  code  method  from  NASA 
Ames  known  as  Panel  Method  Ames  Research  Center  (PMARC),  for  use  in  his  thesis 
research.  The  software  was  accompanied  by  an  additional  routine,  also  developed  by 
NASA  Ames,  known  as  General  Visualization  System  (GVS).  This  second  routine  allows 
a  3-d  graphical  interpretation  of  the  output  from  PMARC,  which  includes  color  display  of 
pressure  distribution,  streamline  generation,  wireframe  modeling  and  solid  surface 
rendering. 

PMARC  enabled  the  designer  to  rapidly  construct  a  complex  geometry  and  analyze 
flows  around  that  geometry  in  the  subsonic  regime.  Since  the  code  allowed  for  assembly 
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of  several  compone  nt  geometries,  it  allowed  for  load  analysis  of  difficult  layouts  that  were 
previously  considered  unworkable  by  most  students  in  the  preliminary  design  phase  (a 
large  fowler  flap  assembly  is  one  example). 

Installation  of  this  code  for  use  by  the  entire  department  did  not  occur  without  some 
additional  effort  beyond  simple  file  transfer  to  the  system  file  server.  The  first  obstacle 
encountered  with  PMARC  was  a  lack  of  adequate  information  for  the  first  time  user 
regarding  the  generation  of  input  and  interpretation  of  output.  To  reduce  the  learning 
requirements  placed  upon  the  user,  an  installation  and  introductory  instruction  manual 
[Ref  S]  was  developed  for  use  by  all  department  members  as  a  part  of  this  thesis  work. 
The  second  hurdle  to  overcome  was  development  of  a  means  to  benchmark  the 
performance  of  the  code. 

1.  Validation  Study  Summary 

This  section  highlights  a  few  of  the  tests  that  were  conducted  to  validate  the 
output  generated  by  PMARC.  Although  one  set  of  validation  results  was  provided  by 
NASA  Ames  with  the  original  documentation,  they  involve  a  rather  complex  wing<body 
combination  that  is  not  well-suited  for  generation  and  duplication  by  the  average  student 
user.  One  example  of  PMARC's  capability  to  produce  an  acceptable  3-d  lift  curve  slope 
using  a  NACA  4  Digit  airfoil  cross-section  is  presented  in  Figure  1 . 

The  comparison  depicted  is  made  from  a  2-d  transformation  method  described  by 
Nicholai  [Ref  3;p.  23]  that  uses  the  turfoil  lift  curve  slope  adjusted  for  overall  geometry. 
Another  comparison  performed  to  verify  the  results  was  on  a  wing  based  using  a 
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NACA  6  Digit  airfoil.  The  comparison  was  made  against  wind  tunnel  data  presented  by 
Bertin  and  Smith  using  an  unswept  trapezoidal  wing  at  low  mach  number  [Ref  6;pp. 
253-257],  Those  results  are  summarized  in  Figure  2. 

Several  validation  studies  were  performed  to  check  the  accuracy  of  output,  in 
addition  to  learning  the  correct  data  input  procedures.  These  studies  demonstrate  a 


Comparison  of  PMARC  and  Transformation  of 
2-D  Section  Data  into  3-D  Piot 


AOA  (degrees) 


3-D  Approximation  PMARC  Solution 

■  « 


Figure  1 

method  for  constructing  wing  shapes  of  various  cross-sections  using  standard  airfoil 
shapes,  in  addition  to  modeling  individually  tailored  designs.  A  complete  summary  of 
these  validation  studies  can  be  found  in  Appendix  A. 

2.  User  Interface  Issues 

Since  PMARC  allows  for  a  wide  range  of  analysis  including  boundary  layer 
investigation,  ground  effect  corrections,  multiple  geometries  and  oscillatory  analysis,  it 
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requires  a  sizable  portion  of  input  data  in  a  FORTRAN  format.  The  code  does  allow  for 
quick  generation  of  any  wing  shape  that  employs  a  NACA  4  Digit  airfoil  by  employing  an 
automatic  panel  generating  scheme.  However  all  other  cross-sections  require  a 
specification  of  the  geometry  in  3-d  space  for  all  dividing  sections;  i.e.,  the  root  and  tip 
cross-section.  To  provide  easier  body  generation,  the  code  allows  several  different  axis 
orientation  systems  to  build  component  geometries.  During  the  development  of  a  user's 
guide,  the  benefit  of  templates  proved  invaluable  in  reducing  data  input  file  construction 
time  and  minimization  of  data  entry  errors.  A  model  template  directory  was  constructed 
on  the  IRIStm  computer  system  that  was  accessible  to  all  users.  This  directory  contained 
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Figure  2 

examples  of  the  construction  of  several  NACA  4  Digit  airfoil  sections  and  one  NACA  6 
Digit  custom  cross-section  design.  With  less  than  a  few  hours  of  instruction  and  practice, 
the  average  user  can  construct  a  NACA  4  Digit  wing  in  approximately  10  minutes.  Given 
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the  coordinates  of  a  custom  cross-section,  most  users  can  construct  an  input  file  using  one 
of  the  templates  in  approximately  20  minutes.  This  short  construction  time  proves  to  be 
very  useful  in  a  time-constrained  design  syllabus. 

For  introductory  design  work,  the  force  and  moment  coefficient  output  of 
PMARC  was  determined  to  be  easiest  to  use  for  most  students  investigating  strength 
requirements  of  a  given  design  proposal.  However,  the  availability  of  pressure  data  offers 
some  interesting  possibilities  for  advanced  study  that  will  be  addressed  later  in  this  report. 
The  lift  and  drag  coefficient  data  provides  information  for  each  individual  panel,  a 
summary  of  all  columns  of  panels  along  the  span,  component  summary  (the  right  wing  for 
example)  and  complete  assembly  summary  in  the  case  of  an  entire  aircraft.  For  this 
project  the  column  summaries  were  used  to  develop  an  approximate  discrete  load  acting  at 
the  centroid  of  each  column  section.  The  column  in  this  output  is  based  upon  the  resultant 
force  and  moment  that  is  obtained  by  grouping  all  of  the  panels  around  a  section  of  the 
airfoil  between  two  neighboring  spanwise  station  locations.  For  instance,  if  the  wing 
panels  divided  the  wing  into  ten  equally  spaced  segments  starting  from  centerline,  the  first 
column  corresponds  to  the  first  segment,  and  all  panels  acting  on  the  upper  and  lower 
surface  of  that  segment.  In  this  analysis,  the  resultant  loads  due  to  lift  and  drag  for  the 
column  have  been  assigned  to  discrete  locations  within  each  column,  i.e.,  the  center  of  the 
column  in  the  spanwise  direction  and  the  center  of  pressure  in  the  chordwise  direction. 
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D.  SAMPLE  PROBLEM  RESULTS 


The  geometry  chosen  for  analysis  is  dq>icted  in  Figure  3,  including  operating 
condition  data.  To  demonstrate  the  procedure  used  in  developing  the  instructional 
modules,  one  aircraft  configuration  was  chosen  to  be  used  as  a  repetitive  example  for  the 
design  sections  of  coursework  encountered  throughout  the  structures  discipline  at  NPS. 
This  allowed  students  to  compare  results  using  different  methods  of  analysis  as  their  skill 
level  and  exposure  to  advanced  methods  increased. 

This  wing  was  chosen  to  support  an  aircraft  operating  at  a  gross  weight  of 
150,000  lbs  subjected  to  a  maximum  allowable  load  factor  of  3.5  g's.  To  demonstrate  the 
methodology  involved,  one  flight  operating  condition  was  used  for  an?«1ysis  and  sizing 
throughout  this  paper,  with  a  discussion  to  follow  on  how  to  apply  these  techniques  to 
multiple  flight  operating  conditions. 

Given  the  allowable  load  factor,  the  condition  chosen  for  analysis  was  slightly  above 
the  comer  velocity.  The  wing  was  assumed  capable  of  delivering  a  maximum  lift 
coefficient  of  1 .35  at  572  ft/sec.  Since  the  panel  code  is  based  upon  potential  flow,  the 

analysis  point  was  moved  out  away  fi’om  c, _ in  an  effort  to  avoid  substantial  bias  due  to 

a  lack  of  non-linear  lift  analysis  capability  in  the  routine.  Using  sea-level  conditions  of 
0.53  mach  and  standard  day  density,  the  required  c^  was  determined  via  equation  (1). 


Cl  = 


2nW 
pV^A^ 


(1) 


The  result  was  a  minimum  of  1 .269  with  the  corresponding  AOA  to  be  determined 
from  the  3-d  lift  curve  versus  Angle  of  Attack  (AOA)  diagram.  To  obtain  a  lift  curve 
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Planform  1 


Cr=\Sft 

C,  =  5ft 
c=  10.83/? 

X  =  0.33 
Aref=  1000/2 
6=  100/ 


Airfoil :  NACA  4418 
^/J=10.0 
««  =  38,378,000 
IMN=  0.5 
Temp  =  60°F 
Altitude  =  Sea  Level 
AOA=2° 


Figure  3 


diagram,  a  PMARC  input  file  for  this  wing  was  created  and  evaluated  at  several  different 
angles  of  attack.  Table  1  is  a  summary  of  those  results  over  various  alphas.  Using  the 
data  of 

Table  1  and  simple  linear  regression  analysis,  the  corresponding  AOA  was  found  to  be 
9.02S  degrees  at  the  desired  operating  condition. 


TABLE  1 


Column 

Q 

-2.0000 

0.2022 

0.0000 

0.3972 

2.0000 

0.5905 

4.0000 

0.7841 

6.0000 

0.9755 

This  angle  of  attack  was  then  used  as  the  input  for  PMARC  to  obtain  the  spanwise 
load.  The  code's  panel  generation  routine  was  employed  for  model  construction.  The 
routine  provided  a  spacing  technique  that  increased  panel  density  near  the  leading  edge, 
trailing  edge  and  wingtips  for  improved  solution  fidelity.  Figure  4  depicts  the  resultant 
column  lift  coefficient  distribution  along  the  semi-span  of  the  right  wing.  In  addition  to 
the  lift  coefficient  output,  PMARC  also  provided  extensive  coefficient  information  in  both 
the  wind  axis  and  body  axis  system.  From  a  structural  standpoint,  the  load  analysis  was 
based  upon  the  normal  and  axial  forces.  Since  these  forces  were  computed  to  act  through 
the  center  of  pressure  on  the  aircraft's  chord  line,  some  modification  of  the  results  was 
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required  to  determine  structural  boundary  conditions  acting  at  the  structure's  centroid. 
Table  2  is  a  summary  of  the  coefficients  that  were  required  along  the  ^an. 


\ 


Spanwise  Lift  Coefficient  Distribution 


V _ y 


Figure  4 

To  simplify  the  following  load  analysis,  the  load  was  assumed  to  be  a  discrete  load 
acting  at  the  center  of  each  column.  To  obtain  the  actual  load,  the  force  coefficients  were 
multiplied  by  the  planform  area  of  each  column  and  the  freestream  dynamic  pressure.  The 
moment  acting  at  each  column  was  obtained  by  multiplying  the  coefficient  by  planform 
area,  dynamic  pressure  and  the  wing  Mean  Aerodynamic  Chord  (MAC,  also  referred  to  as 
cbar).  The  moment  is  the  resultant  torque  due  to  the  normal  force  applied  at  the  column 
center  of  pressure.  The  location  of  moment  computation  was  based  upon  a  user  defined 
reference  line  specified  within  the  input  data  file. 
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TABLE 2 


Ffadmef  p 

o— 

Cs 

Ca 

Cm 

0.0780 

1.2723 

1.2621 

-0.1723 

0.1256 

0.2330 

1.3515 

1.3398 

-0.1877 

0.1934 

0.3820 

1.4066 

1.3932 

-0.2017 

0.2641 

0.5210 

1.4430 

1.4283 

-0.2121 

0.3315 

0.6470 

1.4604 

1.4451 

-0.2171 

0.3903 

0.7580 

1.4492 

1.4343 

-0.2137 

0.4344 

0.8500 

1.3906 

1.3779 

-0.1970 

0.4537 

0.9210 

1.2482 

1.2396 

-0.1619 

0.4330 

0.9690 

0.9951 

0.9933 

-0.1020 

0.3569 

0.9940 

0.7335 

0.7597 

0.0722 

0.2451 

Jj 

_9j 

loj 


III.  STRENGTH  ANALYSIS 


A.  BOUNDARY  CONDITION  DEFINITION 

From  the  discussion  in  Chapter  II  of  the  panel  code  outputs  along  the  semi-qran  of 

the  wing,  it  is  possible  to  obtain  a  reasonable  approximation  of  the  load  pattern.  This 

discussion  will  treat  the  wing  as  a  composite  of  four  bar  elements  as  depicted  in  Figure  S. 

Note  that  these  elements  are  not  constrained  to  be  symmetrical  about  the  X  or  Z  axis. 


Figure  5  also  depicts  the  orientation  of  the  axis  system  that  was  used  throughout  most 
of  the  analysis.  This  sense  of  direction  is  consistent  with  the  output  from  the  panel  code 
where  coefficient  data  was  obtained. 

Since  force  and  moment  computation  requires  knowledge  of  the  reference  area,  the 
planform  area  for  each  spanwise  column  must  be  determined.  The  panel  code  does 
provide  the  location  in  x,y,z  coordinates  for  the  comer  points  of  each  panel  if  instmcted  to 
do  so  within  the  input  file.  From  these  locations  and  simple  trigonometric  reiation^ps. 
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the  pianform  area  of  each  column  can  be  determined  for  use  in  the  calculation  of  applied 
forces  and  moments.  Figure  6  serves  to  describe  the  orientation  of  a  positive  load 
condition.  Since  this  project  was  directed  towards  preliminary,  vice  detail  oigineoing 


design,  the  problem  was  analyzed  at  the  expected  maximum  load  points  near  the  wing 
root.  However,  this  method  can  be  applied  throughout  the  semi-span  to  determine  initial 
minimum  sizing  to  meet  strength  requiremoits. 

At  this  point  in  the  design,  a  common  question  among  new  design  students  pertains  to 
the  choice  of  a  two  spar  system.  Since  this  four  element  structure  is  being  initially 
designed  to  carry  the  bulk  of  normal  stress  due  to  bending,  it  follows  that  a  two  spar 
wingbox  employing  stiffened  wing  skins  will  generally  provide  a  low  cost,  easily 
manufactured,  lightweight  structure.  Another  advantage  as  pointed  out  by 
Niu,  [Ref  7;pp.  247-263]  was  large  interior  volume  for  increased  fuel  capacity.  These 
benefits  seem  well  suited  for  large  commercial  and  military  carriers. 
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This  project  began  by  assuming  the  simplest  possible  configuration  based  upon  a  two 
spar  arrangement.  As  loads  increase  and  size  constraints  decrease  in  a  design  project,  the 
feasibility  of  the  two  spar  structure  begins  to  diminish.  The  designer  needs  to  look  at  his 
results  during  analysis  and  determine  if  they  appear  reasonable.  For  example,  if  the  load 
factor  requirement  is  increased,  the  spar  cap  cross-section  will  also  increase.  Additionally, 
the  wall  thickness  of  the  front  and  rear  spar  wd7S  must  also  be  rising  to  meet  the  increased 
shear  requirements.  As  these  components  grow  in  ^e,  the  designer  needs  to  consider  the 
effect  of  adding  one  or  more  additional  spars  to  begin  reducing  the  structure's  internal 
cross-section  dimensions  (and  ideally  total  wing  wdght). 

Operating  requirements  may  steer  the  designer  in  other  directions,  as  evidenced  in 
current  tactical  military  aircraft  structural  design.  Generally  these  aircraft  demand  high 
strength,  while  also  minimbdng  vertical  wing  depth  to  allow  for  high  speed  operation. 
Most  tactical  aircraft  also  require  multiple  load  paths  for  enhanced  survivability  if 
subjected  to  hostile  fire.  Consider  the  ramifications  of  a  fore  and  aft  spar  structure  that 
has  been  severed  at  the  leading  face  of  the  wingbox  (forward  spar)  under  a  maximum 
loading  condition.  The  unexpected  load  path  may  prove  too  weak  to  carry  the  new  load 
distribution,  resulting  in  catastrophic  failure.  Such  a  design  problem  may  seem 
overwhelming  at  this  point  to  the  beginning  structural  engineer.  Hence,  this  project  has 
been  constrained  to  a  problem  that  readily  allows  solution  via  hand  calculations  within  the 
classroom  v^nthout  reliance  on  finite  element  or  energy  methods. 
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1.  Shear  Diagram  Construction 

Returning  to  the  normal  and  axial  force  coefficients,  the  resultant  loads  acting  at 
each  spanwise  station  of  the  wing  are  depicted  in  Figure  7.  These  loads  provided  the  basis 
for  calculation  of  the  shear  in  the  z  and  x  direction  across  the  face  of  each  section.  Due  to 
the  assumption  that  the  load  acted  at  discrete  points,  the  shear  was  assumed  constant 
across  each  column  elentent.  The  means  of  estimating  shear  along  the  semi-span  was  to 
start  at  the  wing  tip  with  zero  shear  and  sum  the  loads  moving  inboard  to  the  aircraft 
centerline.  Figure  8  represents  the  shear  loads  for  the  two  axes  along  the  semi-^an  of  the 
wing.  For  modeling  considerations,  the  wing  was  assumed  to  be  a  continuous  structure 
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All  shear  values  were  assumed  to  change  at  the  left  face  of  each  column  section  moving 
along  the  y  axis. 

2.  Bending  Moment  Diagram 

The  bending  moment  diagrams  for  the  load  condition  previously  selected  are 
presented  in  Figure  9.  These  diagrams  prove  a  bit  more  challenging  to  construct  due  to 
the  number  of  repetitive  calculations  involved.  Since  the  design  engineer  is  likely  to 
consider  hundreds,  or  even  thousands  of  load  and  geometry  combinations  during  the 
design  sequence,  construction  of  these  diagrams  is  a  prime  candidate  for  automated 
calculation. 

These  diagrams  proved  quicker  to  build  with  the  aid  of  a  spreadsheet  by  starting 
at  the  fixed  end  of  the  wing,  on  centerline.  Figure  10  serves  as  an  example  of  the 
calcuiations  for  the  moment  on  the  left  face  of  column  8  near  the  right  wingtip.  The 
resultant  moment  for  each  face  was  calculated  by  summing  up  the  product  of  all  forces  to 
the  right  of  a  cut  section  and  their  respective  distances  to  the  cut.  These  sums  provided 
the  data  for  the  approximate  moment  at  the  left  face  of  each  column  section. 
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Bending  Moment,  Mxx 


Bending  Moment,  Mzz 


Figure  9 
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Measured  from  Centerline 


Mxx  =  12196(46.05-44.55)  +  5352(48.45-44.55) 
+  1333(49.70-44.55) 

=  46,032  ft-lbf 


Figure  10 
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3.  Torsional  Moment 


Production  of  the  torsional  moment  diagram  depicting  Myy  along  the  semi-span 
requires  the  transfer  of  several  couples  prior  to  the  final  diagram  construction.  PMARC 
provides  a  value  for  the  pitching  moment  coefficient  for  each  column  of  panels  based  upon 
the  planform  area  of  each  column  and  referenced  to  the  Mean  Aerodynamic  Chord  (MAC) 
of  the  overall  wing  structure.  Equation  (2)  serves  to  illustrate  the  relation  between  c^, 

and  the  actual  moment  produced,  M. 

M=^jpl^ArefCCM  (2) 

It  is  important  to  understand  that  c^  and  M  can  be  used  to  specify  moments 

anywhere  on  the  wing  cross-section.  However  the  structural  engineer  is  generally 

interested  in  the  resultant  moments  found  at  the  centroid  of  the  cross-section  under 

investigation.  Anderson  [Ref  7.pp.  lS-30]  provides  an  extensive  discussion  of  the 

integral  approach  to  resolving  the  pressure  force  acting  on  the  surface  of  the  wing  into 

normal  and  axial  components.  These  forces  can  then  be  resolved  into  a  single  pair  of 

forces  acting  at  the  center  of  pressure  on  the  structure.  If  these  forces,  N  and  A  are 

applied  at  the  center  of  pressure  (x^^ ),  there  is  no  induced  pitching  moment  along  the 

chord  line.  Figure  1 1  depicts  the  order  of  tnmsformation  required  to  obtain  a  usable 

relationship. 

Since  the  center  of  pressure  is  unlikely  to  be  located  at  the  cross-section's 
centroid,  a  relationship  must  be  determined  that  can  relate  the  moment  arm  distance 
obtained  between  x^p  and  the  selected  reference.  Returning  to  Anderson's  [Ref  7:p.  30] 
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discussion  of  normal  forces  and  r^ltant  moments,  the  center  of  pressure  is  defined  via 


equation  (3). 


Xcp  = 


Mie 

N 


(3) 


Note  that  this  relationship  only  applies  at  the  Leading  Edge  (LE),  which  is  the  basis  for  the 
transformation  depicted  in  Figure  1 1 .  Anotho-  important  aspect  to  equation  (3)  is  the 
definition  of  positive  orientation  for  the  moment  produced.  Conventional  aerodynamic 
notation  directs  that  any  force  acting  on  the  wing  cross-section  which  causes  a  nose-up 


Figure  11 


rotation  of  the  wing  leading  edge  is  deemed  positive,  (note  that  positive  lift  results  in  a 
negative  M,^.  It  is  also  important  to  realize  that  aerodynamic  sign  covention  for 
moments  and  forces  may  not  coincide  with  the  convention  required  of  structural  analysis. 
It  is  the  design  engineer's  responsibility  to  properly  translate  the  appropriate  sign  when 
moving  among  the  two  disciplines. 

Since  any  moment  can  be  expressed  as  the  sum  of  two  or  more  vector 
components  (the  components  being  moments  themselves),  any  moment  obtained  along  the 
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chord,  where  its  location  and  corresponding  normal  force  is  known,  can  be  used  to  obtain 
the  center  of  pressure.  Since  PMARC  allows  specification  of  the  line  about  which  it 
calculates  c^  (specified  as  RMPX),  and  provides  c^  as  an  output,  it  is  possible  to  locate 
the  center  of  pressure.  Figure  12  highlights  the  relationship  graphically  among  all  relevant 
measurements  for  a  swept  >ving  (referenced  to  the  X=0.0  datum  line).  Equation  (4)  is 
used  to  solve  for  the  location  of  XCP  when  provided  with  the  necessary  input  and  output 
data  from  PMARC. 

XCP  =  mPX-^^^  (4) 

N 

In  equation  (4),  N  is  the  normal  force  obtained  from  the  normal  force  coefficient  c^,  (note 
that  the  output  format  of  the  code  labels  c^  under  the  format  of  Cj^  in  the  body  axis  system 


Figure  12 


and  reserves  c^  as  an  output  of  the  wind  axis  system).  Likewise,  M^^px  moment 
obtained  from  equation  (2)  using  the  output  of  the  code  labeled  c^  (for  the  pitching 
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moment  output  in  the  code  there  is  no  difTerence  betwe«)  in  the  wind  axis  and  the  body 
axis  system). 

The  remaining  translation  required  involves  the  effect  of  the  normal  force  acting 

at  the  center  of  pressure  in  relation  to  the  centroid  of  the  cross-section.  The  resultant 

moment  on  the  centroid  is  readily  determined  as  listed  in  equation  (S). 

mYN  =  NiX^utmi  -  Xcp)  (5) 

Where  is  based  upon  the  location  of  the  centroid,  also  referred  to  as  xbar  in  later 
sections.  Note  that  equation  (S)  specifies  my^  vice  Myy,  this  is  due  to  the  fact  that  the 
chord  line  in  many  airfoils  does  not  pass  through  the  centroid.  To  account  for  this  line  of 
action,  equation  (6)  solves  for  my^,  which  can  then  be  summed  with  my^,  to  yield  the 
desired  resultant  torque  for  an  individual  column  cross-section,  myy.  Once  the  resultant 
moment,  myy,  has  been  obtained  at  each  discrete  column  section,  the  resultant  torsion 
experienced  at  a  cross-section  can  be  determined.  As  previously  depicted  in  the  bending 
moment  discussion,  beginning  at  the  right  wingtip  and  summing  the  column  torsion  values 
yields  the  resultant  torsion,  Myy,  experienced  along  the  span. 

^  W  =  (Zneutml  “  Zc)A 

Due  to  the  format  of  the  output  of  PMARC,  the  torsion  solution  is  obtained  for  each 
column  center,  vice  the  left  face;  recall  that  shear  analysis  was  performed  at  each  column's 
inboard  ftu^e.  Figure  13  is  the  final  result  for  total  toraon  applied  to  the  center  of  each 
column  along  the  ri^t  semi-span  of  the  wing.  Since  Figure  13  was  derived  from  a 
solution  of  equations  (S)  and  (6),  an  interesting  dilenuna  arises  for  the  designer.  Each  of 
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these  equations  assumes  that  the  centroid  is  already  known.  Therefore,  Figure  13  cannot 
actually  be  developed  prior  to  completion  of  the  next  section  involving  solution  of  the 
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Figure  13 

bending  stresses  due  to  N  and  A.  However,  it  is  presented  at  this  time  to  provide  the 
reader  with  a  sense  of  relative  magnitude  in  rdation  to  the  shear  and  bending  diagrams 
previously  developed. 

B.  BENDING  STRESS  ANALYSIS 

The  four  spar  caps  represented  by  the  four  areas  in  Figure  5  will  be  assumed  to 

initially  carry  the  primary  bending  loads  that  give  rise  to  M^x  and  Since  bending 

occurs  in  both  directions,  and  none  of  the  areas  w^e  assumed  to  be  symmetrical,  the 

bending  stress  can  be  found  from  Allen  and  Haider  [Ref  l:pp  166-170],  Assuming  the 

centroid  is  known  and  the  size  of  each  area  is  known,  equation  (7)  serves  to  calculate  the 

normal  stress,  o^,  carried  at  each  spar  cap. 
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Application  of  equation  (7),  in  conjunction  with  an  allowable  working  stress  can  be 


used  to  determine  the  location  of  the  centroid  and  size  of  the  four  spar  caps. 


^  L  /«/»-(/«)* 
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Although  optimization  routines  do  exist  in  industry  today  that  can  quickly  identify  the 
preforred  location  and  distribution  of  the  four  areas  [Ref  8],  observing  the  interaction  of 
the  four  areas  in  relation  to  an  allowable  stress  constraint  proves  very  beneficial  in  the 
learning  environment. 

The  operational  envelope  requirements  set  for  this  exercise  included  a  +3. S  to  -1 .0  g 
load  fiu:tor  boundary.  Analysis  of  the  upper  surface  of  the  wing  reveals  that  at  the 
maximum  load  factor  of  +3.5  g's,  the  structure  will  experience  a  significant  compressive 


load.  Based  upon  information  that  will  be  presented  in  the  buckling  analysis  section,  the 
maximum  allowable  stress  level  for  the  upper  surface  is  the  material  yield  stress  limit. 


which  represents  the  maximum  buckling  stress  for  the  stringers.  However  the  lower  spar 
cap  aUowable  stress  level  for  the  material  chosen  was  set  at  the  ultimate  ten^e  stress. 
Buckling  instability  of  the  lower  surface  did  not  present  the  same  challenge  as  the  upper 
surface  because  of  the  differences  in  design  load  requirements. 

1.  Area  Determination 

As  mentioned  above,  an  iterative  solution  performed  by  the  student  in  design 
proved  very  rewarding  in  understanding  the  demands  placed  upon  various  points  within 
the  cross-section.  Accomplishment  of  this  ta^  was  performed  with  the  aid  of  a 
spreadsheet.  By  initializing  all  four  areas  to  some  arbitrary  amount,  and  then  beginning  a 
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sequence  of  variations,  the  average  student  could  arrive  at  an  acceptable  solution  within 
approximately  ten  minutes  for  a  given  flight  condition.  Appendix  B  provides  a  sample  of 
the  spreadsheet  used  in  this  manner.  Figure  14  is  a  review  of  key  dimensions  employed  by 
the  spar  cap  area  calculation  spreadsheet. 


The  spreadsheet  requires  the  input  of  A1  through  A4,  the  values  of  M^x  and 
and  the  geometry  constraints  of  the  wing  cross-section  at  the  spanwise  location  under 
investigation.  As  the  areas  are  input,  xbar,  zbar,  Izz  and  I2X  are  calculated.  Having 
obtained  these  values,  the  routine  goes  on  to  evaluate  the  stress  at  points  I  through  4, 
using  equation  (7).  While  observing  the  total  area  and  the  stress  levels  for  the  four  points, 
the  designer  brings  points  1  and  3  to  their  respective  allowable  stress  levels,  and  minimizes 
the  total  area  by  varying  areas  at  points  1  through  4.  The  most  productive  method 
discovered  for  minimizing  the  time  involved,  employed  a  rapid  change  of  areas  in  sections 
1  and  3  until  these  sections  are  close  to  meeting  the  appropriate  allowable  stress  levels.  At 
that  point  the  rear  sections  are  varied  until  their  change  be^ns  to  significantly  alter  the 
stress  carried  in  each  respective  forward  member.  At  this  point  it  becomes  important  to 
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start  monitoring  the  total  area  to  determine  which  changes  in  individual  section  areas  leads 
to  a  reduction  in  total  area.  The  process  is  terminated  when  the  solution  with  the  smalls 
total  area  is  found. 

This  problem  was  worked  on  the  basis  of  a  prior  mataial  selection  and  the 
assumption  that  construction  of  stiffeners,  skin  and  spar  caps  would  be  done  with  the  same 
material.  The  method  will  also  work  for  different  material  build-ups  as  long  as  a  nxxlulus 
weighted  moment  of  inertia  approach  is  applied.  The  problem  can  also  be  worked  from 
the  position  of  picking  an  allowable  stress  first  and  then  determining  what  materials  can 
sustain  the  expected  stress  levels. 

2.  Determination  of  Allowable  Stress 

For  those  regions  of  the  structure  that  are  subjected  to  maximum  loads  that  are 
tensile  in  nature,  it  is  of  prime  importance  to  ascertain  what  load  condition  and  stress  level 
should  be  chosen  to  achieve  the  desired  factor  of  safety.  The  following  discussion  was 
applied  to  this  project  in  determining  the  allowable  stress  level  for  maximum  tensile  loads 
for  any  material  where  yield  and  ultimate  strength  test  data  were  available. 

B^in  with  a  comparison  of  two-thirds  of  the  ultimate  tensile  stress  to  the  yidd 
stress  value.  The  smaller  of  these  two  values  will  determine  the  eventual  allowable  stress 
and  limit  load  analysis  conditions,  assuming  a  &ctor  of  safety  of  1 .5.  When  two-thirds  of 
the  ultimate  stress  is  the  lower  value,  analyze  all  loading  problems  at  the  ultimate  load 
(1.5*LimitLoad).  For  this  case  set  g,ii„  .i,i,=g|ihiirmr  However,  if  initially  the  yidd  stress 
value  proved  to  be  the  smaller  value,  the  probimi  should  still  be  analyzed  at  the  ultimate 
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load.  Although,  in  this  situation  the  allowable  stress  level  is  reduced  such  that 
O allowable  =  jG yield  Recall  that  the  limit  load  is  the  maximum  load  used  for  design 

purposes  that  the  structure  is  expected  to  encounter  in  its  normal  operational  routine.  The 

ultimate  load  is  derived  simply  from  the  concept  of  ensuring  the  existence  of  a  margin  of 

safety  in  design  calculations. 

Since  the  upper  si  .iface  of  the  wing  was  subjected  to  large  compression  loads  and 
stringer  buckling  was  the  design  constraint,  the  material  yield  stress  was  chosen  as  the 
allowable  stress.  If,  after  a  buckling  analysis  was  performed,  the  column  length  between 
ribs  was  determined  to  be  unacceptable  (resultant  weight  increase),  one  option  would  be 
to  incrementally  reduce  the  allowable  stress  below  the  yield  value.  This  reduction  in  stress 
carried  will  tend  to  increase  the  weight  of  the  spar  caps,  forcing  a  strength  to  weight  trade 
study  to  compare  spar  cap  increase  in  weight,  to  weight  reduction  due  to  increased  rib 
spacing  (which  results  in  fewer  ribs). 

3.  Spar  Placement  Considerations 

During  the  calculation  of  the  stress  levels  at  each  comer  point  in  the 
cross-section  of  the  root  chord  wingbox,  the  location  of  the  front  and  rear  spars  was 
assumed  to  be  known.  These  parameters  are  not  fixed  at  the  outset  of  the  design  process, 
but  rather  represent  one  additional  consideration  and  potential  optimization  that  will 
require  the  designer's  attention.  Reviewing  Figure  IS  can  provide  a  little  insight  on  the 
considerations  at  hand. 
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Figure  15  depicts  the  location  used  for  the  analysis  in  this  design.  However,  it  is 
important  to  understand  the  range  of  flexibility  involved,  and  associated  implications.  In 
his  presentation  on  layout  techniques,  Niu  [Ref  9;p.  254]  presents  his  range  of  forward 
spar  location  to  be  12-17%  of  the  chord.  He  also  suggests  the  rear  spar  requires 
positioning  at  55-60%  chord  to  accommodate  installation  of  a  30%  aileron  and  flap 
section. 


Thesis  Project 
Spar  Layout 

15%  Chord 


65%  Chord 


Figure  15 

Given  this  range  of  values,  where  does  the  designer  place  the  spar  locations? 
Often  the  aft  spar  is  positioned  based  upon  aerodynamic  considerations  like  attachment 
requirements  to  accomodate  a  flap  or  aileron  system  that  is  necessary  for  control  of  the 
vehicle.  For  this  project  the  rear  wall  of  the  wingbox  was  fixed  at  65%  chord.  Again 
reviewing  Figure  15,  as  the  forward  spar  moves  aft,  the  overall  height  of  the  wingbox  is 
increased  while  the  distance  between  spars  is  reduced.  Locating  the  front  spar  then 
becomes  more  of  a  structural  optimization  problem  that  strives  to  maximize  and 
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to  give  minimum  total  spar  cap  area.  Also,  moving  the  forward  spar  aft,  generally  tends  to 
reduce  overall  wingbox  cross-sectional  area.  This  serves  to  increase  the  shear  flow 
around  the  surface  for  a  given  torsional  loading.  For  the  beginning  designer,  the 
reconunendation  is  to  make  an  initial  comparison  at  the  potential  fore  and  aft  locations  of 
each  spar  to  evaluate  the  magnitude  of  the  differences  found.  For  a  swept  wing  design, 
torsion  nuy  prove  to  be  the  design  driver  as  one  proceeds  out  the  semi-span.  Appendix  C 
offers  further  insight  into  the  options  available  in  a  study  on  various  locations  for  both  the 
front  and  rear  spar. 

4.  Straight  Wing  Approach 

Once  the  initial  effective  area  has  been  determined,  allowing  computation  of  the 
three  moments  of  inertia,  it  is  possible  to  enter  all  required  elements  into  equation  (7)  for 
determination  of  the  effective  stress  at  a  particular  location.  The  only  caution  offered  in 
the  use  of  equation  (7)  involves  maintaining  a  consistent  use  of  orientation.  For  analysis 
of  the  wing  root,  the  evaluation  method  chosen  in  this  attempt  involved  determination  of 
all  forces  and  moments  applied  to  a  positive  face  at  the  desired  analysis  location. 

Figure  16  serves  to  illustrate  the  orientation  chosen.  The  cut  depicted  is  intended  to  occur 
at  an  infinitesimal  distance  from  the  centerline  of  the  wing. 

Since  all  axes  pass  through  the  centroid  for  this  portion  of  the  analysis,  it  is 
important  to  monitor  the  sign  convention  for  both  the  distance  measured  from  the  neutral 
axis  and  the  allowable  stress  levels  used  during  area  determination.  The  sign  orientation 
for  positive  x  and  z  was  based  upon  Figure  16.  The  sign  convention  employed  for  the 
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allowable  stress  assumed  all  compressive  stresses  to  be  negative  and  tensile  stresses  as 
positive.  Appendix  B  also  lists  the  results  of  the  effective  stress  determination  for  several 
operating  conditions. 


Figure  16 


5.  Swept  Wing  Considerations 

One  of  the  key  advantages  in  using  PMARC  for  initial  load  determination  was 
that  it  allowed  analysis  of  virtually  any  geometry  in  wing  design.  So  the  initial  procedures 
for  obtaining  the  loads  applied  to  the  exterior  of  any  structure  operating  subsonically,  will 
always  remain  the  same.  However,  as  Figure  17  illustrates,  application  of  those  loads  will 
dramaticrdly  depend  upon  wing  geometry  chosen.  In  fact,  it  is  a  function  of  both  external 
geometry  (wingsweep  and  dihedral)  as  well  as  internal  geometry  (rib  and  spar  orientation) 
Figure  17  is  taken  from  a  broader  discussion  presented  by  Niu,  where  he  depicts 
the  internal  structural  layout  of  a  wide  range  of  current  commercial  transport  aircraft  that 
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incorporate  swept  wing  designs  [Ref.  9;pp.  278-281],  Since  the  panel  code  output  is 
consistent  with  alignment  shown  in  Figure  17-(a),  one  might  think  that  particular  layout 
would  prove  quicker  to  analyze.  However,  structural  analysis  involves  working  with  an 

I 

axis  of  bending  and  cross-sections  perpendicular  to  it.  The  axis  of  bending  will  coincide 
with  the  line  of  centroids  along  the  span,  which  can  be  taken  as  a  first  approximation  as 

I  Potential  Orientations  I 


(a)  Rib  Aligned  (b)  Rib  placement 

with  A/C  Centerline  perpendicular  to 


Figure  17 

40%  of  the  distance  between  spars.  The  axial  forces  from  PMARC  must  be  resolved  into 
a  component  transverse  to  the  bending  axis  and  along  it  as  shown  in  Figure  18.  The 
transverse  direction  will  be  along  the  rib,  and  the  spanwise  component  will  be  along  the 
stringer.  The  latter  will  add  a  tensile  component  to  the  bending  stresses. 
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To  accomplish  this  task  the  following  methodology  is  proposed; 

•  Resolve  the  loads  obtained  from  PMARC  into  normal  and  axial  componoits 
oriented  with  the  rib  direction. 

•  Modify  equation  (7)  to  account  for  the  axial  load  component.  ' 

•  Update  the  spanwise  position  of  areas  one  through  four  in  stress  calculations. 

Figure  18  hopefully  serves  to  illustrate  transforming  the  applied  load  into  meaningful 
component  requirements.  This  information  provides  the  basis  for  modification  of  the 
analysis  routines  presented  in  Appendbc  B.  At  this  level  of  complexity  in  the  analysis,  the 
amount  of  effort  and  time  consumed  rapidly  begins  to  multiply.  This  should  provide  the 


Figure  18 


reader  with  an  understanding  of  the  requirement  for  more  sophisticated  techniques  for 
analysis  as  the  structure  begins  to  depart  from  conventional  beam  appearance. 
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C.  SHEAR  STRESS  ANALYSIS 

To  allow  for  quick  computation  of  the  expected  shear  on  the  spar  wd>s  and  wing 
covers,  a  lumped  mass  approach,  assuming  constant  shear  flow  between  centers,  was 
employed.  Figure  19  illustrates  the  combination  of  shear  flow  due  to  torsion,  and  shear 
flow  resulting  from  transverse  loadings. 

1.  Shear  Flow  due  to  Torsion 

Since  the  torque  across  the  cross-section  was  solved  in  an  earlier  section  as  the 
resultant  moment  Myy,  the  remaining  solution  requires  minimal  effort.  Equation  (8)  forms 


the  basis  of  a  constant  value  of  shear  flow  that  is  imposed  on  all  faces  between  the  lumped 
mass  centers.  Note  that  equation  (8)  requires  A  to  be  the  cross-section  area  located 
within  the  region  bordered  by  centers  A1  through  A4. 


(8) 


Again,  consistent  sign  convention  is  essential  to  a  correct  solution.  Unlike  bending, 
multiple  sign  errors  throughout  this  analysis  can  result  in  potentially  unrecognizable  final 
solution  errors,  that  in  turn  can  generate  a  structurally  deficient  design.  The  next  chapter 
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will  reveal  the  importance  of  accurate  shear  flow  estimation,  since  serves  as  the  basis 
for  the  buckling  analy»s. 

2.  Shear  Flow  due  to  Applied  Loads 

The  middle  segment  of  Figure  19  also  has  a  straight  forward,  although  longer 
method  of  solution.  For  both  single-cell,  as  the  current  illustration  depicts,  and  multicell 
analysis,  a  set  of  linear  simultaneous  equations  must  be  solved.  To  obtain  the  appropriate 
number  of  equations  in  a  multicell  problem,  the  additional  constraint  of  all  cells  adhering 
to  the  same  angular  rotation  or  deformation  must  also  be  applied.  The  starting  point  in 
determining  %  is  the  relation  between  q.  and  q^,  as  defined  in  equation  (9). 


Alien  and  Haisler  provide  a  discus»on  on  the  derivation  of  equation  (9)  and 
related  multicell  problems  [Ref  1  :pp.  205-228].  Applying  equation  (9)  to  the  middle 
drawing  in  Figure  19,  q,,  q,  and  q,  can  be  related  to  the  shear  flow  q^.  All  that  remains  for 
the  single-cell  problem  is  the  formulation  of  one  last  equation  allowing  solution  of  q,,.  At 
this  p<mt  one  has  the  option  of  analyzing  a  single  face  at  a  given  cross-section,  or  a 
differential  element  containing  a  combination  of  opposing  faces.  Figure  20  has  been 
provided  as  a  simplified  schematic  of  the  options  available.  The  solution  involves  either 
equating  the  moments  produced  by  the  shear  flows  in  the  webs  and  skins  to  the  equivalent 
moments  produced  by  the  resultant  shear  forces  at  that  cross-section;  or  summing  the 
moments  to  zero  for  equilibrium  of  the  shear  flows  and  shear  forces  acting  on  the 
differential  element.  The  results  presented  in  Appendix  B  are  based  upon  analyzing  the 


38 


shear  flow  due  to  concentrated  shear  forces  on  the  same  (positive)  face.  In  this  situation 
the  moments  are  required  to  be  equal. 

The  moments  were  considered  about  A3  since  it  provides  no  torque  effect  from  shear 
flows  qj  and  qj.  Having  solved  for  q^,  equation  (9)  can  again  be  applied,  solving  in 
sequence  shear  flows  q,  through  qj.  These  results  are  algebraically  summed  to  obtain  the 
final  shear  flow  in  each  spar  web  and  wing  cover  (again,  note  the  importance  of  consistent 


orientation).  Figure  21  is  the  final  solution  with  proper  orientation  for  the  operating 
condition  picked  thus  far. 
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Shear  Flow  due  to  Shear  Flow  due  to 

Torsion  Transverse  Shear 

Ih/in  654 


Combined  Shear  Flow 

QRti 


847 


V _ 

Figure  21 


10266 
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D.  VARUTION  IN  OPERATIONAL  CONDITIONS 

The  previous  figures  illustrated  the  shear  flow  developed  due  to  the  operating 

conditions  that  exist  at  point  1  in  Figure  22.  Before  attempting  to  size  the  skin  thickness, 

several  other  operating  conditions  warrant  investigation  to  determine  which  condition 

develops  the  largest  shear  flow  within  the  entire  flight  envelope.  The  four  points 

illustrated  in  Figure  22  encompass  flight  at  high  utd  low  angles  of  attack,  in  addition  to 

maximum  positive  and  negative  load  factors.  In  developing  the  curves  associated  with  the 

maximum  lift,  a  simple  rule  of  thumb  was  employed  to  relate  c,^  of  the  2-d  airfoil  section, 

to  c, _ of  the  3-d  wing  structure.  Obtaining  both  the  positive  and  negative  values  of  c.,. 


for  a  NACA  4418  airfoil  from  Abbot  and  Ooenhoff  [Ref  lO.p.  492],  90%  of  those  values 
were  then  substituted  as  c^^  into  equation  (1)  to  determine  the  corresponding  load 
factors.  The  data  compendium  (DATCOM)  offers  more  extensive  methods  for  estimating 
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the  maximum  lift  of  a  wing  operating  in  the  subsonic  through  hypo-sonic  flight  regimes 
[Ref  11;  Sect  4. 1.3.4]. 

Appoidix  B  provides  ail  stress  levels  ami  shear  flows  for  conditions  I  through  IV. 

For  this  particular  design  problem,  condition  I  was  the  primary  design  constraint  «nce  all 
other  operating  points  at  sea  level  conditions  developed  stress  levels  below  those  of 
condition  I.  If,  in  a  different  design  problem,  one  or  more  flight  conditions  develop  stress 
levels  that  exceed  the  allowable  limits  after  having  azed  the  spar  caps  for  a  given 
condition,  a  logical  sizing  sequence  must  be  employed  to  ensure  the  structure  is  properiy 
sized  to  support  flight  at  all  locations.  The  following  analysis  routine  is  provided  to  aid  in 
that  form  of  design  problem: 

•  Begin  at  condition  I,  note  that  if  PMARC  is  used  condition  I  is  offset  a  small 
amount  fi-om  the  maximum  lift  line  to  avoid  estimation  errors  in  the  nonlinear 
r^ime  of  the  lift  versus  alpha  curve. 

•  Having  sized  the  structure,  move  to  condition  II.  Maintain  all  ^ar  cap  sizes  from 
condition  I  as  the  minimum  values  for  condition  II.  If  the  allowable  stress  is 
exceeded  at  any  point  in  condition  II,  increase  spar  cap  areas  as  required.  Set  these 
new  areas  as  the  minimum  spar  cap  sizes. 

•  Move  to  condition  III;  use  the  same  procedure  as  that  at  condition  II.  Establish  the 
new  minimum  areas  and  move  to  condition  IV.  Having  completed  condition  IV, 
the  structure  is  now  assured  to  carry  the  applied  loads  throughout  the  expected 
flight  envelope. 


E.  SKIN  SIZING 

Based  upon  a  review  of  conditions  I  throu^  IV,  condition  IV  was  found  to  devdop 
the  largest  shear  flow  in  the  upper  skin,  where  q^  was  1238  Ih/in  between  spar  caps.  This 
value  of  shear  flow  was  used  to  determine  the  minimum  ^n  thickness  required  to  carry 
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the  applied  load  along  the  upper  wing  skin.  If  this  value  is  substituted  into  equation  (10) 
and  known,  it  is  possible  to  solve  for  the  minimum  required  skin  thickness  on 

the  basis  of  strength  alone. 


^atlowabtt  “ 


2 

/ 


(10) 


Determining  the  value  of  is  accomplished  in  the  same  manner  as  that  employed  to 

obtain  the  allowable  normal  stress  which  was  calculated  earlier.  To  review  this  method, 
first  consider  the  values  of  Xy^td  and  \Xutt,mtM-  If  the  ItuUimau  proves  to  be  the  lower 

value,  then  set  the  allowable  stress  to  the  ultimate  shear  stress  of  the  material.  However, 

if  the  yield  shear  stress  proved  to  be  the  lower  number,  then  the  allowable  stress  should 
be  set  to  \Xyitu  In  both  of  these  cases,  the  load  used  for  analysis  should  always  be  the 

ultimate  load  using  this  criteria.  For  the  calculations  presented  in  Appendix  B,  material 

properties  for  2024-T4  extruded  aluminum  were  taken  from  Popov  (Ref  12:p.  554], 

however,  a  better  approach  is  to  develop  specific  material  properties  from  the 

MILHANDBOOK-5F  series.  Using  the  information  within  Popov,  the  allowable  shear 

stress  was  taken  to  be  32,000  psi,  which  requires  a  minimum  skin  thickness  of 0.0387  in. 

Having  solved  for  the  skin  thickness  t,  it  is  also  a  responsibility  of  the  designer  to 

investigate  any  other  constraints  placed  upon  the  vehicle  by  outside  agencies  or  the  final 

customer.  Consider  some  of  the  following  requirements  as  set  forth  by  the  Navy.  Table  3 

lists  the  minimum  skin  thickness  requirements  for  a  variety  of  materials  that  may  be 

considered  during  the  design  process  [Ref  13:p.  17].  From  Table  3,  the  minimum  gauge 

requirement  for  alminum  construction  is  0.026  in,  since  this  is  less  than  that  required  due 
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to  the  loads  carried,  the  minimum  thickness  remains  at  0.0387  in.  Had  the  nuuerial 
considered  been  high  strength,  the  minimum  thickness  may  have  been  required  to  increase 
to  comply  with  minimum  gauge  requirements,  this  can  prove  to  be  a  major  factor  in 
considering  a  strength  to  weight  study. 


TABLE  3 


Material 

Mininuim  Gauge  (in) 

Corrosion-reastant  Steel 

0.008 

Alluminum  Alloys  (Exterior) 

0.026 

Titanium  Alloys 

0.016 

Superalloys 

0.015 

Graphite  /  Epoxy  Skins 

0.020 

Boron  /  Epoxy  Skins 

0.020 

One  last  consideration  in  the  manufat^uring  effort  concerns  the  availability  of 
commercial  stock  material.  For  instance,  if  a  roll  of  aluminum  sheet  is  manufactured  in  a 
standard  gauge  thickness  of 0.028  inches,  yet  the  minimum  thickness  required  calls  for 
only  0.026  inches,  what  is  required  of  the  designer?  To  answer  that  question  the  design 
team  would  probably  meet  with  the  material  supplier  to  determine  the  cost  effectiveness  of 
initiating  a  special  production  run,  to  provide  a  limited  amount  of  material  in  the  lighter 
gauge.  The  alternative  to  a  special  production  run  is  acceptance  of  the  associated  excess 
material  (and  weight)  to  achieve  a  lower  manufacturing  cost.  Table  4  is  an  excerpt  from 
Appendix  D  which  lists  some  of  the  dimensions  used  in  production  of  commericial  stock. 
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Note  that  several  standards  of  measurement  still  exist  today,  and  differmt  manu&cturers 
adopt  one  or  more  of  these  standards  for  use  in  their  production  runs. 

In  the  analysis  performed  as  part  of  this  theas,  it  was  assumed  to  be  cheaper  to 
acquire  commercially  available  aluminum  sheet.  Alter  a  review  of  the  load  requirements 
and  the  Navy's  minimum  thickness  requirements  for  aluminum  sheet,  a  minimum  thickness 
of  0.0387  in  is  used  to  enter  Table  4.  With  the  information  available,  it  appears  that  the 
smallest  available  gauge  for  use  in  the  upper  wing  skin  is  on  the  order  of 0.040  in,  which 
corresponds  to  a  number  18  gauge  sheet  in  the  American  StaiKlard.  Using  this  value  as 
the  skin  thickness  t,  and  inserting  it  into  equation  (10)  with  the  orginal  shear  flow,  reveals 
an  effective  working  shear  stress  of  30,948  psi  in  the  upper  wing  cover. 


TABLE 4 


Gauge  No. 

American  Std 

Washburn  Std 

Birmingham  Std 

U.S.  Standard 

18 

0.040 

0.047 

0.049 

0.050 

19 

0.035 

0.041 

0.042 

0.043 

20 

0.032 

0.034 

0.035 

0.037 

21 

0.028 

0.031 

0.032 

0.034 
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IV.  BUCKLING  CONSIDERATIONS 


Having  developed  the  initial  estinute  for  spar  cap  sizes  and  approximated  overall 
shear  flow  in  the  upper  wing  cover,  it  is  now  possible  to  begin  to  estinnate  the  stability 
requirements  of  the  upper  surface.  This  process  takes  on  two  distinct  apects.  Thefir^ 
phase  involves  determining  the  number  and  spacing  of  stiffemrs  to  be  employed.  Once 
that  has  been  accomplished,  the  designer  must  then  choose  the  shape  and  size  of  the 
individual  stringer  elements.  It  is  in  this  portion  of  the  design  that  numy  interactions  occur 
during  calculations  that  open  the  door  to  a  wide  range  of  trade-offs  that  the  designer  must 
review.  This  diversity  of  possibilities  in  the  deagn  is  perhaps  one  of  the  most  challenging, 
yet  rewarding  aspects  of  the  process  for  the  apprentice  designer. 

A.  STIFFENER  SPACING 

Recall  that  during  the  bending  analysis,  the  allowable  normal  stress  level  for  the  upper 
suTflu:e  was  set  at  the  material  yield  stress.  One  should  also  observe  that  the  upper  cover 
deals  with  two  distinct  forms  of  primary  buckling  aiudysis,  colunm  buckling  of  the  stiffener 
and  plate  buckling  for  the  skin  itself  The  plate  buckling  problem  becomes  the  starting 
point,  where  the  operating  stress  in  the  skin  is  kept  below  the  buckling  stress  of  the  skin 
when  treated  as  a  flat  plate.  The  methodology  presented  here  is  based  upon  an  attempt  to 
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preclude  buckling  throughout  the  wing  structure  in  all  flight  conditions  operated  at  or 
below  ultinute  load  conditions. 

The  working  shear  stress  now  becomes  the  defining  buckling  shear  stress  for  the 
wing  skin.  The  dimensions  of  the  plate  shown  in  Figure  23  will  be  selected  to  prevent 
shear  buckling.  The  equation  for  shear  buckling  of  flat  plates  is  given  in  equation  (11) 
[Ref  2;p.  172].  The  desired  output  of  equation  (11)  is  the  ratio  of  plate  thickness  to 
width. 


Xcr  ~ 


(11) 


j^E 
112(1 

since  the  critical  buckling  stress  is  known,  plus  the  ^n  thickness,  t,  it  is  possible  to  solve 


-V>)JUJ 


for  plate  width,  b.  The  plate  width  is  the  distance  between  stringer  attachment  lines. 
Referring  to  Figure  23,  if  the  stringers  are  attached  with  a  single  row  of  closely  spaced 
rivets,  then  b  would  represent  the  distance  between  these  rivet  lines,  which  are  depicted  by 
dashed  lines  along  the  span. 


Figure  23 
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Before  one  solves  for  b  explicitly  in  equation  (II),  several  key  constants  have  to  be 

obtained.  The  first  is  the  modulus  of  Elasticity,  E.  The  results  presented  in  Appendices  B 

and  C  were  based  upon  2024  aluminum  possessing  a  uniform  modulus  of  elasticity  of 
10.6x10*^  Ib/sq-in  for  aluminum  plate.  The  constant  v  is  Poisson's  ratio,  and  was  taken  to 
be  0.33.  The  final  constant,  K„  is  known  as  the  shear  buckling  constant  for  use  with  the 

Euler  buckling  equation.  The  value  of  K,  reflects  the  type  of  support  the  plate  has  at  its 

boundaries.  It  ranges  from  5.35  for  simply  supported  edges  to  8.98  for  four  rigidly  held 

edges  in  plates  that  are  rectangular  in  shape  [Ref  2;p.  172].  Since  ideal  boundary 

conditions  are  seldom  the  rule  in  actual  structures,  review  of  several  design  handbooks  can 

provide  a  range  of  values  obtained  from  experience  to  use  for  K,.  Niu  provides  several 

illustrations  of  design  charts  for  K,  that  depend  upon  how  each  side  of  the  plate  is 

supported,  the  ratio  of  the  plate  length  to  width  and  the  curvature  of  the  surface  of  the 

plate  [Ref  9:pp.  138-140]. 

For  simplicity  of  calculations,  the  value  chosen  in  this  analysis  was  based  upon  the 
discussion  of  Sechler  and  Dunn  in  their  development  of  shear  buckling  for  a  flat  plate. 

They  provide  a  figure  that  theoretically  estimates  this  constant  for  two  different  end 
supports.  One  curve  is  based  upon  simply  supported  sides,  while  the  second  employs 
clamped  edge  support.  Since  the  actual  aircraft  skin  panel  support  structure  falls 
somewhere  between  these  two  locations,  an  average  was  taken  for  K,,  assuming  the  length 
of  the  plate  (which  will  eventually  determine  the  rib  spacing)  would  be  at  least  5  times 
longer  than  the  width  of  the  plate.  For  these  parameters  the  value  of  K,  was  found  to  be 
7.17.  [Ref  2;pp.  164-173] 
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Determination  of  the  total  number  of  stringers  required  along  the  upper  cover 
necessitates  the  knowledge  of  the  plate  width  between  stringer  attachment  lines.  Since  the 
locations  of  the  leading  edge  spar  and  trailing  edge  spar  have  already  been  fixed,  the 
distance  between  them  is  also  known.  If  the  stringers  are  placed  at  equal  distances 
between  these  two  spars,  equation  (12)  can  be  used  to  solve  for  total  number  of  required 

stringers  for  the  upper  wing  cover  between  the  two  spars.  Note  that  the  number  of 

Toni  Width  _| 
b 

individual  plate  elements  exceeds  the  number  of  stiffeners  by  one  if  the  fore  and  aft  edges 
are  assumed  to  be  bounded  by  the  spar  caps  themselves.  Solving  for  flight  condition  IV, 
b=l  .42  in  and  n=62  total  stringers  at  the  wing  root. 

B.  STIFFENER  SIZING 

This  portion  of  the  buckling  analysis  offers  many  possible  variations  to  arrive  at  an 
acceptable  solution.  Rather  than  provide  the  student  with  an  optimal  solution,  it  was  again 
determined  to  be  more  beneficial  for  the  student  to  explore  the  effects  of  variations  in  a 
wide  range  of  parameters.  The  first  obvious  question  to  arise  in  sizing  of  a  stiffener  stems 
fi'om  the  geometry  of  the  stringer  itself  Niu  is  quick  to  point  out  the  importance  of  the 
design  engineer  understanding  more  than  optimum  strength  or  minimum  strength.  In  his 
description  of  the  Y  stringer,  he  indicates  that  it  is  a  very  strong  support  that  can  be  easily 
attached  during  manufacture.  However  this  same  device  is  expensive  and  provides  an 
inaccessible  region  on  the  inside  skin  surface  that  may  become  susceptible  to  corrosion 
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damage  during  the  life  of  the  aircraft.  This  region  also  precludes  a  visual  inspection 
without  complete  removal  of  the  stringer.  Another  concern  of  the  designer  is  the  inter&ce 
of  numufticturing  processes  involved  in  physically  attaching  a  stringer.  Questions  to  be 
asked  include,  "Does  this  design  hinder  the  shop's  ability  to  rivet  or  attach  jig  and 
alignment  hardware?"  Awareness  of  these  factors  at  the  earliest  point  in  the  process  can 
produce  enormous  savings  over  the  life  of  a  program.  [Ref  9;pp.  141-142] 

Keeping  these  ideas  in  mind,  the  Z  section  was  chosen  to  demonstrate  one  method  of 
stringer  sizing.  The  Z  geometry  does  not  place  enormous  computational  requirements  on 
the  designer,  is  readily  attached  and  allows  for  visual  inspection. 

1.  Initial  Size  Estimate 

In  an  attempt  to  incorporate  previous  design  experience  and  reduce  the  time  for 
solution,  two  important  ratios  were  taken  from  Niu  as  minimum  weight  rules  of  thumb  and 
applied  to  the  stringer  cross-section  build-up.  The  first,  listed  in  equation  (13),  ratios  the 
total  cross-sectional  area  of  the  stringer  to  the  cross-section  of  the  plate  located  between 
two  stringers. 

^,  =  1.5{A/)  (13) 

The  second  relation  in  equation  (14)  is  taken  specifically  for  use  with  the  Z  cross-section 
geometry  and  relates  the  skin  thickness,  t,  to  the  thickness  of  the  vertical  web  of  the 
stringer,  t,. 

/,  =  1.05/  (14) 

Returning  to  one  last  manufacturing  consideration,  actual  construction  of  the 
stringer.  One  more  cost  study  that  may  prove  beneficial  early  in  the  process,  is  the  choice 
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of  stringer  build-up.  Possible  options  include  commercial  stock,  limited  production 
extrusions  or  bending  operations  to  produce  the  desired  section.  If  it  is  determined  to  be 
cost  effective  to  bend  a  flat  sheet  into  the  appropriate  Z  section,  one  should  observe  that 
the  value  taken  for  t,  now  becomes  constant  for  all  three  segments  of  the  Z.  Turning  to 
Figure  24,  it  depicts  one  possible  labeling  of  the  required  dimensions  used  in  this  analysis 
to  arrive  at  the  final  geometry. 


Figure  24 


Also  observe  the  segment  defined  by  2We  in  Figure  24.  The  value  of  2We  is 
taken  from  the  concept  of  an  efifective  width  of  skin  acting  in  support  of  the  attached 
stringer  to  prevent  its  compression  buckling.  It  should  not  be  confused  with  the  width  of 
the  plate  between  stringers,  which  is  greater  than  the  effective  width.  In  Figure  25,  the 
general  pattern  of  stress  variation  can  be  seen  as  one  moves  in  a  chordwise  manner  along 
the  wing  cover.  As  expected,  the  stress  level  rises  in  the  vicinity  of  the  stringers,  which 
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must  resist  compression  buckling.  The  stress  carried  then  falls  off  to  its  low  point  midway 
between  stringers  where  only  the  skin  acts  to  carry  the  resultant  load. 


Figure  25 


The  effective  width  relates  the  stress  carried  at  the  supported  edge,  labeled  Om, 
to  the  normal  buckling  stress  of  the  plate  section,  referred  to  as  the  critical  stress,  Ocr. 
Since  the  upper  wing  cover  structure  was  initially  sized  to  ensure  the  maximum  nornud 
stress  level  remained  below  the  material  yield  stress,  the  value  of  Om  is  taken  to  be  o>«u. 
In  the  case  of 2024  aluminum,  the  yield  stress  was  assumed  to  be  44,000  psi.  Using 


anotho*  form  of  the  Euler  buckling  equation,  the  critical  buckling  stress  for  a  flat  plate  can 
be  used  to  determine  (Jcr,  as  described  in  equation  (15).  Again,  the  values  of 


will  vary  depending  upon  the  type  of  edge  ^ipports  the  flat  plate  encounters.  The  values 


of  K,  range  from  4.0  for  the  simply  supported  case,  to  7.5  for  all  sides  rigidly  clamped. 


Again,  taking  an  average  of  these  values,  K,=5.75  was  used  for  the  solution  of  equation 


(15)  in  determing  the  critical  stress  of  the  plate  sections  between  stringer  attadiments. 


52 


Substituting  this  value  into  equation  (IS)  results  in  a  buckling  stress  of  24, 181  psi  for 
flight  condition  IV  in  the  upper  skin  plates  along  the  wing  root.  It  should  be  noted  that 
the  largest  value  the  critical  stress  can  become  is  the  material's  normal  yield  ^ress. 

[Ref  2:pp.  164-169] 

Returning  to  the  relation  of  the  critical  stress  of  a  plate  section  and  the  stress 
carried  at  its  edge  supports,  equation  (16)  defines  this  relationship  in  terms  of  the  effective 
width  of  the  skin  attached  to  one  stringer.  Sechler  and  Dunn  provide  an  extensive 
discussion  on  the  development  of  equation  (16)  and  determination  of  the  period  of 
oscillation  associated  with  the  loading  pattern  depicted  in  Figure  25  [Ref  2;pp.  203-234]. 

Using  these  inputs,  the  effective  width  is  readily  obtained  as  0.74  in  or  1 .48  in  for 
2We  as  depicted  in  Figure  24.  This  information  is  now  used  in  another  Euler  buckling 
problem  that  can  determine  the  ratio  of  column  required  length  to  radius  of  gyration  for 
the  section  shown. 
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2.  Uaifomi  Tkkkncss  Approach 


If  the  stringer  thickness  is  assumed  to  be  uniform  throughout  (as  in  the  case  of  a 

section  produced  by  bending  up  a  flat  sheet),  a  starting  point  in  estimating  the  required 

»ze  is  quickly  obtained.  Since  the  overall  area  of  the  stringer  was  determined  via  equation 

(13)  and  the  thickness  in  equation  (14),  it  becomes  possible  to  rdate  the  total  length  of  the 

Z  stringer  sub-sections  in  the  form  of  equation  (17). 

/l,  =  (/l+/2  +  /3)/,  (17) 

Turning  to  prior  design  experience,  the  relative  dimensions  of  the  legs  that  produce  a 

structure  for  optimum  strength  are  specified  in  equations  (18)  and  (19). 

A3=0.3(/2)  (18) 

/3=0.7(/l)  (19) 

Once  the  length  of  each  segment  is  known,  it  becomes  possible  to  locate  the 
center  of  gravity  for  the  entire  stringer  cross-section  that  is  depicted  in  Figure  26.  This 


Figure  26 

composite  stringer  includes  the  skin  section  of  2We  and  all  three  stringer  segments. 
Following  determination  of  the  centroid,  the  next  requirement  is  to  calculate  the  moment 
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of  inertia,  1„  for  the  composite  cross-section.  It  is  the  combination  of  effective  skin  width 
and  the  attached  stringer  that  will  act  in  unison  as  an  individual  column  subjected  to 
compressive  buckling  along  the  span  (between  the  wing  ribs).  Recalling  that  the  moment 
of  inertia  is  related  to  the  area  of  the  cross-section,  as  illustrated  in  equation  (20),  it 


becomes  possible  to  solve  for  the  radius  of  gyration,  r. 

s-Le- 

A  total 

All  that  remains  for  this  first  approximation  is  determination  of  the  required 


(20) 


colunrm  length  or  rib  spacing.  For  buckling  beyond  the  proportional  limit.  The  tangent 
modulus  is  inserted  into  the  Euler  column  curve,  and  it  is  possible  to  develop  equation 
(21),  which  can  be  used  to  solve  for  the  column  length,  L. 


Ocr  .  _ f  Ocr  ^ 

-7~rv 


( 


jl] 

Jet) 


(21) 


The  constant  m  in  equation  (21)  represents  the  inverse  of  the  strain  hardening  exponent 
for  the  given  material.  In  this  analysis  Of 2024-T4  aluminum,  m  was  S.  Similarly,  K  also 
represents  a  material  property,  it  is  known  as  the  strength  coefficient.  It  should  be  noted 
that  K  is  generally  provided  in  units  of  ksi;  this  requires  consistent  matching  of  units 
among  the  values  substituted  into  equation  (21).  Throughout  this  analysis  a  value  of 
K=1 17  ksi  was  used. 

3.  Stringer  Cross-Section  Stability 

Although  the  results  obtained  above  were  based  upon  approximate  size 
relationships  obtained  from  previous  design  experience,  it  is  still  important  to  verify  the 
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integrity  of  the  stringer  itself  subjected  to  ultimate  loading.  The  easiest  means  of  analysis 
is  accomplished  by  treating  each  individual  segment  of  the  Z  stringer  as  a  plate  demrat 
subjected  to  compressive  loading.  Again  returning  to  the  Euler  buckling  rdation,  this  time 


in  the  form  of  compressive  plate  loading,  results  in  equation  (22). 

.  K„k^E  (i,\‘ 


<Tcr,= 


The  value  of  Kc,  varies  depending  upon  the  location  of  the  indidual  strings  dement. 

Referring  to  Figure  24,  note  that  the  upper  flange  and  lower  flange  pieces  (described  as 

parts  1  and  3),  are  assumed  to  act  as  a  simply  supported  plate  on  three  sides,  with  the 

fourth  side  free.  Using  the  plate  normal  stress  buckling  results  provided  by  Sechler  and 

Dunn  [Ref  2;p.  168],  reveals  K,=K]=K).SO.  Similarly,  the  web  section  (part  2),  is 

considered  to  be  a  simply  supported  on  all  four  sides;  thus  K2=4.00. 

Now  substituting  these  values  into  equation  (22)  allows  calculation  of  the 

expected  buckling  stress  for  each  individual  plate  element  composing  the  stringer.  For 
each  element,  the  result  is  then  compared  to  the  yield  stress.  The  lower  of  Ojhm  and  Oer  is 

then  referred  to  as  ciauowabu,  which  is  used  to  find  the  crushing  stress,  Cee,  for  the  entire 
stringer.  To  determine  this  crushing  strength  begin  by  finding  the  allowable  effective  load 
in  each  element  given  by  equation  (23). 

(23) 

Once  the  allowable  load  in  each  element  is  determined,  the  crushing  stress  is  determined 


via  equation  (24). 


Oce  = 


(Pith+PA 


(Ai+A2+Ay)  '  ' 

The  importance  of  this  check  is  to  ensure  that  the  crushing  strength  of  the  stringer  is  in 
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&ct  greater  than  the  material  yield  stress.  If  the  value  comes  up  short,  then  the  designer  is 
forced  to  change  the  values  of  A,  through  A,  until  that  goal  is  met.  Recall  that  if  the 
manufacturing  intent  is  to  construct  the  stringer  from  bending  a  single  sheet  of  material, 
this  places  a  restriction  on  area  changes  due  to  length  changes  only.  However,  if  custom 
extrusions  are  being  considered,  then  both  length  and  thickness  of  the  individual 
components  are  available  for  customizing  to  meet  design  strength  requirements. 
Considering  just  these  few  variables  that  are  open  to  design  selection  at  only  one  operating 
condition,  begins  to  portray  the  expansive  nature  of  optimization  involved  with  an 
aircraft's  structural  design  life  as  thousands  of  variables  surface  for  use  in  trade  studies. 
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V.  STRENGTH  TO  WEIGHT 


This  chapter  offers  the  beguining  designer  some  ideas  and  options  of  trade  studies 
that  can  be  used  to  optimize  the  previous  deagn.  These  ideas  are  but  a  few  of  the  many 
available.  The  distinction  between  this  academic  project,  and  actual  manufiKture  of  an 
aircraft,  arises  from  the  need  to  devdop  and  understand  principles  on  a  small  scale  prior  to 
their  application  in  the  6eld.  The  topics  will  be  presented  in  the  order  of  bending  stress 
considerations,  followed  by  instability  constraints.  However,  the  order  of  these  studies 
also  occur  in  a  chronological  sequence  that  supports  the  current  sequence  of  teaching 
within  the  department. 

A.  BENDING  STRENGTH  TRADE  STUDY 
1.  Material  Selection 

Since  each  material  has  a  unique  set  of  ultimate  and  tensile  load  properties,  the 
value  of  GtOowabit  varies  by  a  large  amount  over  different  types  of  conventional  metals. 

Referring  to  equation  (7)  and  inserting  different  values  for  the  allowable  stress  will 

generate  significant  variations  in  the  resultant  areas  that  exist  in  the  four  spar  ciq>s.  For 

each  case  considered,  every  material  has  a  unique  density  associated  with  that  material. 

For  ease  of  comparison,  consider  a  wing  section  at  any  given  location  along  the  span  and 

assume  it  to  be  of  unit  length  in  the  y  direction.  For  instance,  in  this  problem  consido-  a 

chordwise  segment  that  includes  the  wing  from  y=K>.00  out  to  y=l  .00. 
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Next  assume  the  cross-section  areas  of  the  four  spar  caps,  (A1  -  A4),  remain 
uniform  throughout  the  entire  section  (y=0.00  :  y=l  .00).  Consider  the  case  of 
2024  aluminum  alloy,  assumed  to  weigh  0. 100  lb  /  cu-in.  For  the  problem  solved  in 
Chapter  III,  the  total  weight  of  the  spar  caps  for  a  unit  depth  at  the  wing  root  is  indicated 
by  equation  (25). 

Wcaps  =  (A  1  +A2+A3  +^4)p  •  ^gravi^  (25) 

Note  that  in  equation  (25),  the  density  is  assumed  to  be  in  units  of  mass  /  cu-ft.  When  the 

density  is  specified  in  that  fashion,  an  adjustment  for  dimensions  must  be  included  to  arrive 

at  the  correct  units  for  total  weight  per  unit  span.  For  the  density  value  listed  above,  the 

acceleration  due  to  gravity  has  already  been  included.  Table  5  is  a  quick  summary  of 

representative  values  for  a  few  common  engineering  materials. 


TABLES 


Material 

Weight  /  Volume 
(Ib/in') 

Ultimate  Tensile 
Strength  (ksi) 

Yield  Strength 
(ksi) 

2024-T4  Aluminum 

0.100 

60 

44 

606 1 -T  6  Aluminum 

0.098 

42 

36 

5Cr-Mo-V  Steel 

0.281 

240 

200 

17-7PH  Stainless 

0.276 

170 

140 

Titanium 

0.162 

115 

110 

Inconnel  718 

0.297 

180 

150 

One  should  also  notice  the  lack  of  two  additional  columns  in  Table  5,  these 
include  the  cost  per  pound  to  acquire  the  material  and  a  similar  reference  of  cost  per 
finished  foot.  Consider  the  case  of  titanium,  for  certain  applications,  the  tooling  required 
to  shape  a  particular  component  may  prove  unacceptable  in  the  overall  cost  constraint  of 
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the  project  for  a  limited  production  run.  The  average  student  is  unlikely  to  possess  such 
detailed  cost  data,  but  should  understand  its  impact  when  progressing  to  an  industrial 
design  team  or  contract  review  team. 

2.  Forward  Spar  Location 

Earlier  in  chapter  III,  it  was  pointed  out  that  the  front  spar  location  would 
generally  occur  between  12  -  17%  of  the  chord.  Since  there  is  an  available  range  to 
anchor  the  forward  spar,  this  offers  another  opportunity  to  evaluate  the  effect  of  changing 
the  spar's  location.  A  series  of  spar  locations  might  be  investigated  for  a  given  operating 
condition,  aircraft  configuration  and  material  selection  to  determine  the  optimum  location 
of  the  forward  spar.  Again  the  optimum  point  will  be  identified  by  the  lowest  total 
cross-sectional  area  required  of  the  spar  caps.  An  example  of  this  type  of  investigation  has 
been  provided  in  Appendix  C. 

B.  INSTABILITY  IMPLICATIONS 

The  shear  flow  in  the  skin,  and  thus  minimum  required  skin  thickness,  will  change 
similar  to  the  change  in  spar  cap  size  given  a  variation  in  material  choices.  However,  the 
shear  flow  problem  offers  an  additional  consideration  not  associated  with  bending.  The 
earlier  analysis  finished  with  solving  for  the  skin  thickness  that  sustained  the  ultimate 
expected  load  and  complied  with  Navy  design  specifications.  Now,  the  designer  needs  to 
consider  the  effects  of  increasing  the  skin  thickness  beyond  the  required  minimum.  As 
skin  thickness  rises,  so  does  the  load  carried  by  the  stiffeners.  This  increase  tends  to 
reduce  the  column  length  in  response  to  increased  buckling  potential.  However,  as  the 
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buckling  column  lengths  are  reduced  (implying  reduced  rib-spacing,  thus  more  ribs 
required  along  the  span),  the  total  number  of  stiffeners  required  is  also  reduced 

Aside  from  changing  the  skin  thickness  directly,  rib  spacing  can  also  be  affected  by 
the  value  specified  for  the  upper  surface  allowable  stress  design  level.  Additionally,  for 
any  given  thickness  chosen,  if  the  stringer  crushing  strength  proves  to  be  unsatisfactory, 
the  areas  of  the  individual  stringer  elements  can  be  varied.  All  of  these  changes  impact  the 
ratio  of  skin  thickness  to  rib  spacing.  As  both  of  these  quantities  change,  so  does  the 
overall  weight  of  the  wing.  For  educational  purposes  it  is  recommended  that  at  least  one 
trade  be  conducted  by  the  students  and  investigated  over  a  reasonable  range  of  parameters 
to  introduce  the  concepts  and  gain  an  understanding  of  the  extent  of  possibilities  involved 
over  the  entire  flight  envelope  of  an  aircraft. 

C.  ADVANCED  MATERIAL  CONSmERATIONS 

As  a  student's  understanding  and  exposure  to  structural  considerations  grows,  one 

more  consideration  that  weighs  heavily  in  all  aviation  designs  remains  for  review.  The 

design  engineer  must  consider  the  effects  of  cyclic  loading.  This  arena  may  force  the 

complete  exclusion  of  various  materials  for  failing  to  meet  minimum  required  life  cycles. 

As  a  student  completes  his  or  her  study  in  the  area  of  fracture  and  fatigue,  he  or  she  can 

then  analyze  effectively  the  likelihood  of  a  given  design  satisfying  the  Navy's  expected  life 

cycle  requirements.  This  stage  of  the  design  could  prove  to  be  the  most  challenging  in 

attempting  to  improve  the  strength  to  weight  ratio  of  the  structure  while  subjected  to  a 
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wide  range  of  loads  over  a  long  lifetime.  Structural  fatigue  has  been  a  central  issue  in  all 
Naval  Aviation  programs  for  the  past  30  years. 
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VI.  SUMMARY 


A.  FUTURE  AREAS  OF  STUDY 
1.  Aerodynamic  Investigations 

During  the  validation  and  testing  of  PMARC,  all  of  the  work  was  completed 
without  the  use  of  the  boundary  layer  routines  that  are  built  into  the  analysis  portion.  This 
was  done  to  prevent  overwhelming  a  student  in  his  or  her  first  exposure  to  such  a  complex 
code.  The  boundary  layer  routine  may  prove  however  to  offer  an  interesting  and  effective 

means  of  determining  the  approximate  c, based  upon  a  reversal  in  the  boundary  layer 

shear  stresses.  A  separate  research  project  in  itself  would  encompass  a  validation  of  the 
boundary  layer  routine  and  eventual  incorporation  into  the  user's  guide  that  was  developed 
as  part  of  this  thesis  project. 

Since  most  high  lift  devices  are  employed  well  below  Mach  0.5,  PMARC  remains 
a  good  solution  for  those  unusual  design  problems.  However,  the  supersonic,  transonic 
and  high  subsonic  regimes  cannot  be  addressed  currently  through  the  employment  of 
PMARC.  These  regions  comprise  a  major  portion  of  the  operational  envelopes  of  today's 
tactical  jet  aircraft.  To  support  the  structures  design  of  wing  sections  operating  at  high 
speeds,  with  all  high  lift  devices  stowed,  a  simplified  CFD  routine  could  be  developed  for 
student  use.  Such  an  undertaking  would  require  earlier  introduction  of  CFD  methods  into 
the  current  syllabus,  a  grid  generation  system  that  is  readily  adaptable  to  geometric 
iterations,  and  analysis  routines  that  can  respond  to  large  variations  in  freestream 
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conditions.  Having  met  all  of  those  conditions,  the  output  needs  to  be  transformed  into  a 
usable  format  for  structural  analysis,  in  a  fashion  similar  to  that  presented  in  this  thesis,  or 
like  the  pressure  distribution  analysis  described  in  the  next  section.  Given  the  magnitude 
of  that  challenge,  the  distinction  between  educational  and  industrial  goals  needs  to  be 
addressed  at  length  prior  to  undertaking  such  a  project. 

2.  Advanced  structural  Analysis 

Several  areas  of  improvement  were  identified  throughout  this  project  and 
presented  for  consideration.  The  first  point  introduces  the  incorporation  of  Finite  Element 
Analysis  (FEA)  into  the  advanced  aircraft  structures  course.  This  is  the  first  introduction 
to  the  subject  for  most  students  in  the  aerospace  discipline  at  NFS.  A  logical  building 
block  on  the  original  design  project  would  be  verification  and  comparison  of  results  using 
the  layout,  sizing  and  material  choices  from  the  AA2021  class  project. 

If  the  initial  class  project  were  to  be  modified  to  include  design  and  analysis  of  a 
swept  wing,  as  depicted  in  Figure  27,  the  wing  could  be  analyzed  over  the  period  of  two 
quarters.  The  first  quarter  would  be  done  using  the  methods  presented  here  with  an 
adjustment  for  loading  based  upon  the  leading  edge  sweep.  The  root  section  would 
require  the  use  of  FEA.  The  root  section  analysis  would  also  necessitate  the  incorporation 
of  the  results  of  the  straight  wing  approach  applied  at  the  boundary  of  the  two  regions. 
This  project  probably  represents  both  the  most  extensive  of  follow-on  suggestions,  and  the 
most  beneficial  to  future  design  students. 
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If  fewer  fundamental  concepts  are  presented  in  future  instruction  to 
accommodate  increased  FEA  skills,  there  are  two  investigation  areas  to  aid  in  expanded 
use  of  FEA.  The  first  is  oriented  towards  the  reduction  of  time  consumed  to  generate  a 
finite  element  model.  Since  PMARC  is  already  used  to  determine  a  wirefi'ame  model, 
incorporate  the  means  to  port  that  information  into  a  FEA  routine  as  a  part  of  an  input 
file.  PMARC  is  written  in  FORTRAN  77;  thus,  only  a  minor  modification  should  be 
required  within  the  source  code  to  generate  an  additional  output  data  file.  This  project 
would  only  be  available  to  students  of  U.S.  citizenship  however,  due  to  restrictions  on 
accessing  the  source  code. 


Figure  27 


The  second  project  entails  the  boundary  conditions  that  are  solved  under  FEA. 
As  mentioned  earlier,  PMARC  can  generate  output  on  the  basis  of  pressure  distribution 
(Cp)  or  force  and  moment  coefficients.  Since  most  FEA  routines  can  solve  an  external 
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pressure  problem,  a  modification  to  the  source  code  to  provide  yet  another  output  data 
file,  could  be  used  to  input  the  boundary  conditions  to  a  particular  FEA  application.  This 
project  would  parallel  the  previous  one  in  development,  with  only  slight  differences 
resulting  from  the  requirement  of  converting  Cp  to  actual  pressure. 

B.  CONCLUSION 

The  computer  code  developed  to  produce  the  outputs  displayed  in  Appendbc  B 
represents  a  significant  step  forward  for  the  aircraft  design  course  offered  to  NFS 
students.  It  offers  a  means  to  reduce  required  computational  time  without  sacrificing  the 
importance  of  careful  consideration  of  the  impact  of  engineering  decisions  in  establishing  a 
baseline  design.  Hopefully  the  methods  presented  in  the  body  of  this  thesis  in  conjunction 
with  the  sample  calculations  listed  in  Appendix  C  offer  insight  to  the  beginning  design 
student  that  can  build  one's  confidence  level  in  that  first  design  undertaking. 

The  intent  for  the  software  development  on  a  spreadsheet  format  was  to  reinforce  the 
student's  understanding  in  a  guided  environment.  As  each  portion  of  the  series  of 
spreadsheets  becomes  relevant  in  a  student's  classroom  discussion,  exact  portions  are 
worked  out  in  the  classroom.  Ideally  the  project  will  support  about  four  hours  of  in-class 
instruction  (for  all  core  structures  courses),  to  allow  the  student  the  opportunity  to  derive 
the  underlying  equation  used  in  every  calculation  found  throughout  the  worksheets. 
Additionally,  the  classroom  exercise  should  introduce  the  concept  of  working  in 
engineering  teams  to  achieve  a  final  result  in  a  reasonable  amount  of  time.  By  constraining 
the  project  to  an  in-class  basis,  the  instructor  can  ensure  participation  by  all  students  and 


provide  valuable  insight  on  stumbling  blocks  as  they  arise  among  individuals.  This 
intermediate  approach  to  design,  should  serve  students  extremely  well  as  they  prepare  for 
eventual  placement  in  the  most  demanding  acquisition  jobs  the  Navy  can  offer,  both  today 
and  on  the  horizon  ahead. 
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APPENDIX  A 


A.  WAKE  INITULIZATION 

The  panel  code  used  for  this  analysis,  PMARC,  allows  for  a  wide  range  of 
computations.  One  of  its  most  complex  features  includes  the  analy«s  of  flow  past  an 
oscillating  airfoil.  To  accomplish  this  task,  the  code  requires  the  assembly  and  analy«s  of 
a  wake.  Unfortunately,  that  type  of  computation  generally  demands  extensive  CPU  time 
on  most  computers.  Recognizing  that  handicap,  the  creator  of  PMARC  wanted  to  allow 
experienced  users  of  the  code  several  means  by  which  they  could  reduce  the  overall  run 
time  during  program  execution,  depending  upon  their  individual  requirements.  As  a  result 
the  user  has  the  option  to: 

*  Specify  the  number  of  timesteps  over  wluch  the  wake  is  calculated.  This  assumes 
analysis  b^ns  at  time  ^K),  where  the  body  under  investigation  is  about  to  enter 
the  freestream  flow. 

*  Specify  the  length  of  the  time  incremoit  used  for  completion  of  one  time  step. 

•  Allow  PMARC  to  initialize  the  trailing  wake  and  determine  its  devdopment  in 
<lirection  and  nuignitude. 

•  Override  the  code  and  allow  the  user  to  spedfy  the  wake's  initial  aze  and  direction 
upon  startup.  This  option  is  designed  to  reduce  overall  run  time  if  the  user  is 
simply  trying  to  analyze  a  complex  flow  geometry  at  a  steady-state  condition. 

Because  of  the  availability  of  these  two  features,  and  the  intent  of  developing  several 
templates  that  would  satisfy  a  wide  range  of  student  applications,  a  determination  of  the 
optimum  settings  for  the  student  design  requirement  had  to  be  accomplished.  To  review 
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the  effects  of  varying  the  number  of  timesteps,  a  test  matrix  was  conducted  on  a  consistent 
wing  geometry  to  observe  when  the  outputs  converged  within  a  reasonable  tolerance  to  a 
steady-state  value.  In  reviewing  Figure  28,  the  steady-state  solution  was  found  to  occur 
to  within  four  decimal  place  agreement  following  10  timesteps.  This  test  was  performed 
using  the  geometry  described  as  Planform  3  in  Reference  S.  The  operating  conditions  for 
the  test  were  the  defaults  listed  for  Planform  3  and  AOA=0.00  degrees. 


Variation  in  Number  of  Timesteps. 
Convergence  to  a  Steady-State  Soiution 


Location  on  tha  Right  Wing  (%  Sami-Span) 


1  TinUp  STmalv*  10  Timiiapi  ISTinmlips 
-•  ♦  •  -B- 

Th«  Mikl  line  rapieMnIs  the  anal  coiweiyenoe  folowing  10  Mmeetepe 


Figure  28 

Turning  to  the  question  of  which  initialization  method  to  employ,  the  user 
specification  of  wake  attachment,  or  allowing  PMARC  to  initialize  the  wake,  prompted 
the  development  of  a  second  test  scenario.  The  first  concern  in  this  test  questioned  the 
existence  of  reasonable  agreement  between  a  user  defined  wake,  and  PMARC  initiation  of 
the  wake.  Since  the  wake's  departing  orientation  is  generally  quite  complex,  it  will  not 


69 


immediately  adopt  the  direction  of  the  freestream,  which  can  lead  to  errors  in  the  analyas. 
Following  that  result  came  the  question  of  the  existence  of  any  time  savings  afforded  to 
the  typical  NFS  student  using  the  various  initiation  options.  The  final  review  would  then 
focus  on  the  accuracy  of  the  results.  Figure  29  is  useful  to  demonstrate  the  first  point. 
Using  the  same  airfoil  fi’om  the  earlier  test,  the  two  methods  of  wake  generation  were 
compared  at  AOA=0.00  and  AOA=8.00  degrees. 


Comparison  of  Wake  Generation  Methods 


BAHAMA 


0*0  + 
PMARCMfiiirtoirrtesdupon^T^ 

User  Sp«cMad  VMi0  bnad  upon  4  TlmMtafM 
foHcMNng  lh«  NtW  slart^ip  calcuM^ 


Fig^e  29 

Reasonable  agreement  was  obtained  between  the  two  methods,  therefore  a  second 
test  was  conducted  to  determine  the  time  savings  involved  by  specifying  the  wake's  initial 
direction  and  magnitude.  Figure  30  is  a  summary  of  those  results.  At  first  inspection  of 
Figures  29  and  30,  one  might  assume  a  distinct  advantage  exists  in  user  specification  of 
the  wake  initial  positioning.  However,  it  should  also  be  noted  that  the  code  rendered 
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good  convergence  for  the  problem  after  only  10  timesteps.  At  that  iteration  level,  the  time 
savings  offered  by  the  user  specification  mode  is  not  considerable.  An  additional  point 
involves  how  much  time  is  required  to  actually  develop  the  required  information  that  is 
specified  in  the  input  data  file  to  accomplish  user  initialization.  From  that  aspect,  coupled 
with  the  intent  to  develop  templates  to  be  available  for  general  use,  the  determination  was 
made  to  allow  PMARC  to  specify  the  initial  wake  orientation  and  magnitude.  The  penalty 
of  additional  run  time  does  not  outweigh  the  advantage  of  significant  potential  errors  that 
may  be  induced  due  to  user  specification  of  the  wake. 


B.  ANALYTICAL  VALIDATION 

Having  established  the  desired  input  format,  the  next  step  in  the  process  involved 
checking  the  output  against  other  methods  of  calculation.  The  verification  test  first 
attempted  to  match  the  value  of  total  lift  for  a  wing  against  a  prediction  based  upon  a 
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conversion  technique  described  in  Nicholai.  Knowing  the  characteristics  of  a 

two-dimensional  airfoil  (infinite  span),  Nicholai  offers  a  means  to  convert  those 

characteristics  into  a  finite,  three-dimensional  wing.  Equations  (26)  and  (27)  were 

employed  to  accomplish  this  conversion  of  data.  In  these  equations  the  constant  C 1  was 
assumed  to  be  0.6,  while  the  value  for  ^  was  taken  as  0.866.  Note  that  P  is  a  correction 

for  compressibility  which  is  not  considered  within  PMARC.  Equation  (27)  is  actually  the 

first  to  be  solved  since  it  determines  the  lift  curve  slope.  [Ref  3;pp.  2-8;  2-9] 


Cl  =  [C£«]j^(a  -Oiouft)  +  Cl(a  -olouji) 
2nAR _ 


Cu  = 


2+ 

The  results  of  this  comparison  are  depicted  in  Figure  3 1 .  Since  the  approximation 
accounts  for  compressibility,  it  shows  a  slightly  larger  value  of  lift  for  any  given  AO  A  at 
the  operating  condition  of  O.S  Mach. 

Having  attained  a  good  correlation  between  total  lift  and  AOA,  the  next  task  was  to 
look  at  the  lift  distribution  the  code  predicted  along  the  span  at  a  given  AOA.  Planform  3 
was  again  chosen  as  the  baseline  configuration  for  the  test.  The  analytic  model  chosen  for 
the  reference  was  based  upon  a  modification  of  Schrenk's  Lift  Approximation  as  described 
by  Raymer.  The  process  begins  by  first  solving  for  the  average  between  a  predicted  chord 


fl  +  tan!A' 

'  PM 


(26) 

(27) 
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length  determined  via  equation  (28)  and  the  actual  (or  trapezoidal)  chord  length  at  various 
stations  along  the  span.  The  station  points  chosen  were  increments  of  10%  of  the 
semi-span. 


The  value  of  S  is  obtained  from  equation  (29),  note  all  other  parameters  are  taken  from 
the  characteristics  ofPlanform  3. 

5  =  |C(1+X)  (29) 

The  trapezoidal  chord  length,  C(y),  at  a  particular  station  location,  y,  can  be  found  from 
equation  (30). 

C(y)  =  c|l-^(1-X)j  (30) 


Figure  31 

Once  the  elliptical  chord  length  and  trapezoidal  chord  length  have  been  determined, 
Schrenk's  Approximation  states  that  the  distribution  of  total  lift  along  the  span  is 
proportional  to  the  average  of  these  two  chord  lengths.  If  this  average  chord  length  at  a 
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given  station  is  plotted  against  the  station  value,  the  total  area  under  the  resultant  curve 
must  sum  to  be  equal  to  the  total  lift.  [Ref.  4;pp.  342-343] 

To  accommodate  a  comparison  between  PMARC  and  Shrenk's  ^proximation,  a 
somewhat  different  curve  was  developed  using  the  same  principles.  Once  the  average 
chord  values  were  determined  for  the  approximation,  they  were  normalized  by  dividing 
each  individual  value  by  the  root  chord  length.  Since  this  approximation  was  based  upon 
lift,  vice  lift  coefficient,  the  output  of  PMARC  also  required  modification.  Taking  the 
local  coefficient  of  lift  for  each  spanwise  station,  that  value  would  have  to  be  multiplied  by 
the  dynamic  pressure  (q)  and  the  area  of  that  spanwise  station.  Having  obtained  the  lift  at 
each  station,  all  values  were  then  normalized  by  dividing  the  lift  at  the  first  output  location 
(7.8%  of  the  semi-span)  of  PMARC  into  all  station  lift  values. 

Since  the  lift  distributions  for  both  methods  were  based  upon  a  location  described  by 
percentage  of  semi-span,  it  was  then  possible  to  overlay  the  two  results  for  accurate 
comparison.  These  results  are  depicted  in  Figure  32.  For  reference.  Table  6  provides  a 
brief  synopsis  of  the  characteristics  as^iociated  with  Planform  3. 


TABLE  6 


1  Planform  3  Characteristics  | 

-b 

36.7  ft 

c, 

20.0  ft 

0.1445 

420  ft* 

^le 

25  d^ 
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Normalized  Load  Carriage 


Location  on  the  Right  Wing  (%  Semi-Span) 


Figure  32 
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C.  COMPUTER-BASED  VERIFICATION 


Having  attained  reasonable  results  from  generalized  approximations,  the  next  step 
entailed  comparing  the  output  of  PMARC  with  an  independent  panel  code.  The  choice 
for  this  comparison  was  a  panel  code  named  SUB,  which  was  developed  by  NASA 
Langley.  The  test  involved  nearly  identical  geometry  and  freestream  conditions. 

However,  two  small  differences  did  exist  between  the  two  codes.  First,  SUB  did  not 
permit  easy  identification  of  the  wing's  cross-section,  as  a  result  the  airfoil  considered  in 
SUB  was  a  flat  plate.  The  cross-section  used  for  PMARC  was  a  NAC A  00 12  airfoil.  The 
second  difference  was  due  to  the  selection  of  independent  variable  used  by  each  code  in 
their  analysis  routines.  In  SUB,  the  user  specifies  the  desired  total  output  c,^,  whereas 
PMARC  evaluates  the  input  on  the  basis  of  operating  angle  of  attack.  Figure  33  depicts 
the  actual  wing  chosen  for  this  portion  of  validation  testing. 


Figure  33 


The  results  of  those  tests  are  illustrated  in  Figure  34.  The  runs  were  compared  along  the 
right  wing  for  two  different  choices  of  AOA.  Even  with  the  subtle  differences  in  the 
codes,  the  results  demonstrated  good  correlation. 
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Figure  34 


D.  EXPERIMENTAL  DATA  COMPARISON 

The  last  validation  study  included  a  check  against  actual  wind  tunnel  results  for  a 

wing  section  that  was  built  with  a  NACA  6  Digit  airfoil  cross-section.  Data  for  the  study 

were  obtained  from  Benin  and  Smith  in  their  discussion  of  panel  code  methods 

[Ref  6:pp.  253-256],  The  test  cross-section  was  a  NACA  65-210  forming  a  15  ft  span 

wing  similar  in  shape  to  Pianform  1  used  throughout  the  body  of  this  thesis.  Comparison 

of  the  wind  tunnel  results  and  PMARC  output  is  displayed  in  Figure  35.  The  effects  of 

non-linearity  due  to  viscous  effects  can  be  observed  at  the  higher  values  of  AOA.  With 

the  conclusion  of  this  test,  PMARC  was  determined  to  be  a  good  introductory  tool  for  use 

by  design  students  at  NPS,  that  could  also  be  readily  validated  by  available  references. 
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APPENDIX  B 


This  appendix  provides  the  output  of  the  spreadsheet  routines  developed  to  support 
design  students  faced  with  the  task  of  repetitive  calculations  based  upon  the  principles 
presented  in  this  thesis.  The  results  provided  are  for  flight  conditions  I,  II,  III  and  IV, 
which  were  previously  described  within  chapter  III  of  this  thesis.  These  values  are  based 
upon  the  analysis  of  an  aircraft  wing  geometry  represented  by  Planform  I  throughout  the 
thesis.  The  routines  used  to  determine  normal  and  axial  loadings  allow  the  inclusion  of 
planform  drag  through  the  specification  of  c^.  The  values  obtained  at  all  four  flight 
conditions  used  the  airfoil  section  data  at  the  appropriate  velocity  and  Reynold's  number 
from  Abbot  and  von  Doenhoff  [Ref  10;p.  493].  Anderson  provides  an  example  of  using 
the  section  drag  characteristics  to  solve  for  the  expected  lift  and  drag  associated  with  a 
finite  wing  [Ref  14:pp.  221>222]. 

Each  p«^e  represents  a  complete  output  from  individual  worksheets  that  have  been 
linked  together  in  a  3-D  spreadsheet  format.  The  code  was  developed  within 
LOTUS®  1-2-3  release  4.0  for  WindowSp^.  The  code  runs  without  delay  on  a  386-based 
machine  or  better  that  employs  a  math  co-processor  and  at  least  8  Mb  of  RAM.  For  those 
users  with  a  color  VGA  or  better  monitor,  the  input  cells  have  been  shaded  to  further 
reduce  the  time  required  to  conduct  a  preliminary  design  study.  Those  inputs  dealing  with 
a  wing's  physical  geometry  are  indicated  in  light  blue.  Yellow  cells  pertain  to  material 
properties,  while  purple  is  used  for  those  items  that  should  always  be  updated  with  any 
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change  to  the  operating  conditions.  Due  to  the  reproduction  limitations  of  this  thesis, 
within  Appendix  B,  aU  input  cells  are  indicated  by  a  grey  background. 


Panel  Code  Planform  Area 
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APPENDIX  C 


This  appendix  was  developed  to  demonstrate  the  process  of  design  variation  for 
optimization,  more  commonly  known  as  a  trade-study.  This  appendix  will  consider  two 
specific  cases  that  support  study  within  the  current  core  structures  curriculum  at  NPS. 
The  first  study  is  intended  for  use  within  the  introductory  structures  course,  while  the 
second  study  is  aimed  at  the  advanced  structures  student. 

A.  LOCATING  THE  FORWARD  SPAR 

Earlier  in  Chapter  III,  it  was  presented  that  the  forward  spar  location  will  generally 

fall  between  the  12%  and  17%  chord  position,  assuming  a  conventional  two  spar 
configuration.  Since  this  range  of  decisions  is  available  to  the  designer,  he  or  ^e  must 
determine  the  final  location.  Since  the  average  introductory  level  engineer  lacks 
considerable  design  experience,  a  reasonable  and  quantitative  process  must  be  devised  to 
lead  one  to  that  final  decision. 

One  recommendation  is  to  begin  first  by  reviewing  what  parameters  have  been  used 
as  inputs  in  determing  the  strength  of  the  wing  box.  Considering  the  spar  cap  area  sizing 
spreadsheet,  one  observes  that  several  inputs  are  affected  by  changing  the  location  of  the 
forward  spar.  These  inputs  include  XI,  X2,  Z1  and  Zc.  Through  the  use  of  a  scale 
drawing  of  the  wing's  cross-section,  it  is  possible  to  locate  the  inputs  required  for  a  wide 
range  of  ^)ar  locations. 
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However,  before  entering  a  long  trial  and  error  process  of  arbitrarily  entering  a  series 
of  spar  locations  into  the  spreadsheet,  a  short  review  of  anticipated  outputs  can  prove  very 


beneficial  in  reducing  computation  time  and  increasing  one's  understanding  of  the 
principles  at  work.  One  or  more  specific  outputs  should  be  identified  as  being  uniquely 
related  to  specific  inputs.  For  this  problem,  the  first  obvious  output  to  consider  is  total 
area  of  the  r  oar  caps.  Since  maximum  strength  in  comparison  to  structural  weight  is 
sought,  the  smaller  the  total  area  at  a  fixed  allowable  stress,  the  greater  the  overall 
strength-to-weight  ratio.  However,  moving  the  structure's  centroid  may  also  affect  the 
shear  flow  due  to  torsion  around  the  body.  As  shear  flow  increases,  so  does  skin 
thickness,  which  serves  to  increase  total  weight.  Table  7  is  a  summary  of  data  for  moving 
the  forward  spar  within  the  range  of  12%-17%  of  chord. 

TABLE? 


Spar  Location 

Total  Area  (in^ 

Shear  Flow,  (Ib/in) 

12% 

152 

998 

13% 

150 

1,007 

14% 

148 

1,015 

15% 

147 

965 

16% 

142 

1,018 

17% 

142 

960 

As  observed  in  Table  7,  moving  ihe  spar  aft  to  17%  reduces  the  total  area  without  a 
major  fluctuation  in  shear  flow  in  the  upper  wing  cover.  Returning  to  the  discussion  on 
anticipation  of  results.  Consider  the  NACA  4418  airfoil  used  in  the  wing  cross  section. 
For  the  12%-17%  chord  range,  as  one  moves  aft  a  sizeable  increase  in  the  wingbox  height 


occurs.  This  implies  that  zbar  is  also  increasing  in  the  calculation  of  the  bending  stresses. 
Yet,  xbar  is  decreasing  at  the  same  time,  so  why  the  inconsistancy  in  results?  The  answer 
lies  in  the  sensitivity  to  change  of  the  two  parameters.  The  chord  change  from  12%  to 
13%  produces  a  relative  increase  in  ^ar  compared  to  the  height  of  the  forward  spar  that  is 
much  greater  than  the  reduction  in  xbar  compared  to  the  overall  length  of  the  wingbox. 

This  revelation  should  lead  the  designer  to  two  additional  considerations.  The  first 
issue  that  should  come  to  mind  is  the  quesion  of  consistency.  Does  this  trend  in  improved 
structural  performance  hold  for  moving  the  forward  spar  aft  for  all  airfoils?  Hopefully,  the 
beginning  student  would  respond  no,  since  the  rate  of  vertical  change  to  horizontal  change 
is  a  function  unique  to  every  airfoil  design.  Second,  what  other  assumptions  were  made  in 
the  development  of  Table  7?  Recall,  the  rear  spar  was  fixed  at  6S%  chord.  That  fact 
should  warrant  at  least  on  investigation  of  the  points  in  Table  7  at  a  different  rear  spar 
location.  A  good  place  to  start  might  be  the  forward  limit  of  travel  for  the  rear  spar, 
which  is  generally  taken  to  be  50  %  of  the  chord. 

B.  RIB  SPACING  CONSIDERATIONS 

In  Chapter  IV  of  this  thesis,  a  method  was  developed  based  upon  previous  design 

experience  to  optimize  a  stringer  size  and  location  for  material  strength.  When  the  issue 
of  total  weight  enters  that  process,  often  very  different  design  decisions  will  result  that  will 
change  one's  recommendations  rather  drastically.  Recall  that  as  the  rib  spacing  is 
increased,  the  stringer  cross-sectional  area  must  also  increase  to  preclude  column  buckling 
of  the  stiffener.  Also  an  increase  in  the  wing  skin  thickness  reduces  the  effective  working 
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shear  flow  carried  in  the  skin.  As  the  shear  flow  drops,  the  distance  between  stringers 
increases.  This  distance  increase  will  imply  a  reduction  in  the  total  number  of  stringers, 
and  hence  lower  total  stringer  weight. 

These  points  should  lead  the  designer  to  consider  the  effect  of  increasing  skin 
thickness  on  total  weight  of  the  wing.  Does  an  optimum  condition  exist  that  will  carry  the 
applied  loads  for  lowest  total  wing  weight?  Additionally,  what  input  parameters  must  be 
considered  in  such  a  study?  These  types  of  questions  are  fundamental  for  efficient  and 
effective  work  in  the  field  of  engineering  design. 

To  determine  the  input  parameters,  consider  those  variables  in  the  current  structure 
that  add  weight  to  the  wing.  The  components  considered  thus  far  include  the  spar  caps, 
the  wing  skin  and  the  ribs.  Knowing  the  density  of  the  material  and  the  dimensions  of 
each  structure,  it  is  possible  to  estimate  the  weight  of  each  of  these  components.  In  the 
case  of  the  spar  caps,  the  cross-section  area  times  the  length  of  the  semi-span  yields  the 
volume.  Yet  this  volume  does  not  change  when  the  skin  thickness  is  varied,  therefore  it 
does  not  enter  into  this  particular  trade  study. 

For  a  quick  example  of  trade  study  calculations  that  are  relevant,  the  following 
assumptions  are  made  for  analysis; 

•  Only  the  wing  skin  on  the  top  of  the  wing  between  the  forward  and  rear  spar  will  be 
considered  for  total  wing  weight  estimation. 

•  The  initial  study  will  consider  a  wing  of  constant  chord,  equal  to  the  root  chord 
dimension  (a  rectangular  wing). 

•  The  only  wing  skin  associated  variable  affecting  total  weight  will  be  skin  thickness. 
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Note  that  to  properly  determine  the  location  of  the  rib  spacing,  the  actual  wing  chord 


at  various  spanwise  stations  must  be  accounted  for.  Turning  to  the  only  remaining 
variable,  some  estimation  process  must  be  derived  for  the  weight  of  the  ribs.  For  the 
results  depicted  in  Figure  36,  the  following  steps  were  used  to  calculate  individual  rib 
weight; 


•  Determine  the  cross-sectional  area  bounded  by  the  four  spar  caps  to  use  as  the 
external  boundary  of  the  rib. 

•  Assume  only  70%  of  that  area  contains  sheet  metal,  that  allows  for  cutouts  and 
weight  reduction. 

•  The  rib  weight  is  the  product  of  the  reduced  cross-sectional  area,  material  density 
and  specified  rib  skin  thickness.  For  illustration  purposes,  three  different  skin 
thickness  were  used,  beginning  with  the  minimum  gauge  value.  Actual 
determination  of  the  required  rib  skin  thickness  is  generally  a  detail  design 
procedure. 

Figure  36  is  the  result  of  the  spreadsheet  that  follows  at  the  end  of  this  appendix.  The 
calculations  are  based  upon  the  same  equations  developed  in  chapter  IV  of  this  thesis. 
However,  one  distinction  is  made  for  this  method  that  differs  from  the  previous 
development.  Recall  earlier,  the  rib  spacing  (L),  was  determined  by  first  specifying  the 
stringer  cross-section  geometry,  and  then  solving  for  and  the  radius  of  gyration.  When 

L  is  specified  first,  the  designer  must  then  solve  for  the  radius  of  gyration  and  required 
nranwnt  of  inertia.  The  geometry  must  then  be  developed  to  satisfy  the  resultant  moment 
of  inertia  and  radius  of  gyration.  This  method  appears  to  require  a  bit  more  iteration  on 
the  part  of  the  designer,  but  does  lead  to  an  overall  reduced  weight.  Hopefully,  these 
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studies  offer  some  insight  into  the  nature  of  optimization  that  serves  to  build  a  skill  that 
can  be  applied  to  a  wide  range  of  engineering  decisions. 


Total  Upper  Wing  Weight  vs  Rib  Spacing 


Assum«8  uniform  chord  length 
Only  accounts  for  skin  waigM  on  upper  wing  surface 


Figure  36 


Rib  Spacing  Trade  Study 
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